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1.0  INTRODUCTION 


In  missile  preliminary  design  it  Is  necessary  to  quickly  and  economically 
estimate  the  aerodynamics  of  a  wide  variety  of  missile  configuration  designs. 
Since  the  ultimate  shape  and  aerodynamic  performance  are  so  dependent  upon  the 
subsystems  utilized,  such  as  payload  size,  propulsion  system  selection  and 
launch  mechanism,  the  designer  must  be  capable  of  predicting  a  wide  variety  of 
configurations  accurately.  The  fundamental  purpose  of  Missile  Datcom  is  to 
provide  an  aerodynamic  design  tool  which  has  the  predictive  accuracy  suitable 
for  preliminary  design,  and  the  capability  for  the  user  to  easily  substitute 
methods  to  fit  specific  applications. 

The  Missile  Datcom  computer  code  calculates  the  static  stability  and  con¬ 
trol  characteristics  using  the  methods  discussed  in  the  Missile  Datcom  Hand¬ 
book.  Although  the  companion  Handbook  is  not  required  to  use  the  computer 
code,  it  serves  as  documentation  for  the  incorporated  methods. 

This  manual  consists  of  two  volumes.  This  Volume  is  the  program  User's 
Guide.  Instructions  for  incorporating  the  program  on  the  user's  computer  sys¬ 
tem  are  detailed  in  Volume  2.  The  program  has  been  written  to  closely  conform 
to  the  American  National  Standards  Institute  (ANSI)  Fortran  IV  language,  and 
is  written  for  the  Control  Data  Corporation  (CDC)  computers.  The  code  is  com¬ 
patible  with  the  Fortran  V  language  and  is  available  in  a  VAX  version.  It  can 
be  easily  converted  to  other  computer  systems. 

Questions  regarding  the  code  should  be  directed  to  the  United  States  Air 
Force  Wright  Aeronautical  Laboratories,  AFWAL/FIGC,  Wright-Patterson  AFB,  Ohio 
45433,  Telephone  (513)  255-8485. 
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2.0  PROGRAM  CAPABILITIES 

The  computer  code  is  capable  of  addressing  a  wide  variety  of  conventional 
missile  designs.  For  the  purposes  of  this  document,  a  conventional  missile  is 
one  which  is  comprised  of  the  following: 

0  An  axisymmetric  or  elliptically-shaped  body 

0  One  to  four  fin  sets  located  along  the  body  between  the  nose  and 
base.  Each  fin  set  can  be  comprised  of  one  to  eight  identical  panels 
attached  around  the  body  at  a  common  longitudinal  position 

0  An  airbreathing  propulsion  system. 

To  minimize  the  quantity  of  input  data  required,  commonly  used  values  for 
many  inputs  are  assumed  as  defaults.  However,  all  program  defaults  can  be 
over-ridden  by  the  user  in  order  to  more  accurately  model  the  configuration  of 
interest.  .  / 

/ 

The  following  paragraphs  detail  the  configurations  that  can  be  analyzed. 
Later  paragraphs  describe  the  range  of  aerodynamic  coefficients  that  can  be 
predicted.  Finally,  the  program  constraints  are  discussed. 

2.1  Addressable  Configurations 

The  following  configurations  can  be  analyzed: 

a)  Circular  or  elliptically-shaped  cross  section  bodies,  with  or  without 
airbreathing  inlets 

b)  Fin  alone  (1  to  8  panels  attached  at  the  root) 

c)  Body  and  up  to  4  finsets  (1  to  8  panels  in  each  finset) 

d)  The  body  and  finset  configurations  with  deflected  fins 
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Certain  restrictions  exist  due  to  method  limitations  and  are  summarized  in  the 
following  paragraphs. 


2.1.1  Axi symmetric  or  Elliptical  Bodies  -  Methodology  is  incorporated  that 
permits  analysis  of  the  configuration  components  summarized  in  Table  1.  Due 
to  the  types  of  methods  selected  restrictions  also  apply  to  the  manner  in 
which  thei.=*  components  are  joined  to  form  a  complete  configuration: 

Subsonic/transonic  speeds  -  The  aerodynamic  methods  assume  that  the  body 
is,  as  a  minimum,  composed  of  a  nose-cylinder  combination.  The  afterbody 
(boattail  or  flare)  is  optional,  but  if  used,  it  must  be  attached  to  a 
cylindrical  center  body  whose  length  is  at  least  four  body  diameters;  this 
restriction  minimizes  nose  flow  field  coupling  over  the  afterbody.  If  an 
afterbody  is  specified  it  must  not  be  cylindrical,  e.g.,  the  base  diameter 
must  be  different  than  the  centerbody  diameter.  Table  2  summarizes  the 
other  restrictions  on  the  configurations. 

Supersonic  speeds  -  The  aerodynamic  methods  used  are  not  restricted  to 
nose-cylinder  combinations.  Any  arbitrary  radii  distribution  can  be  de¬ 
fined  since  theoretical  techniques  are  employed  at  Mach  numbers  above 
1.4.  Care  should  be  taken  to  avoid  Introducing  unexpected  corners  into 
the  contour.  If  the  contour  has  any  concaved  regions  the  marching  may 
fail  due  to  shock  impingement  on  the  body  as  it  starts  to  curve  out. 

2.1.2  Panels  -  The  program  will  accept  inputs  to  describe  most  airfoil  sec¬ 
tions  or  planforms.  Certain  assumptions  and  limitations  are  made  and  summa¬ 
rized  in  the  following  paragraphs. 

2. 1.2.1  Airfoil  Section  -  The  program  will  accept  virtually  any  symmetrical 
airfoil  section  or  NASA  subsonic  cambered  section.  The  airfoil  section  can  be 
defined  using  a  NACA  designation  or  by  supplying  the  coordinates  of  the  sec¬ 
tion.  Circular  arc,  hexagonal,  or  diamond  shaped  sections  can  also  be  speci¬ 
fied.  A  symmetric  hexagonal  cross-section  is  the  default;  its  shape  is  com¬ 
puted  using  the  planform  inputs.  Hence,  explicit  definition  of  the  airfoil 
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section  Is  optional.  Although  cambered  airfoil  sections  can  be  Input,  their 
use  In  the  code  Is  currently  limited  to  subsonic  applications, 

2.1. 2. 2  Planform  -  Each  set  of  fins  may  be  comprised  of  up  to  eight  separate 
panels.  It  Is  assumed  that  each  panel  Is  geometrically  Identical.  Although 
planforms  may  be  described  by  up  to  10  separate  pieces  or  sections,  an  equiva¬ 
lent  straight-tapered  panel  1s  computed  and  used  at  all  speeds.  There  is  no 
capability  to  specify  a  panel  with  outboard  dihedral. 

2.1.3  Airbreathinq  Inlets  -  Both  axisymmetric  and  two-dimensional  airbreath¬ 
ing  Inlets  can  be  defined  on  the  geometry.  They  can  be  In  combinations  of  1, 
2,  3  or  4  Identical  Inlets  and  can  be  positioned  around  the  body  at  arbitrary 
angles.  Due  to  strong  component  Interference  effects,  inletted  vehicles  can 
only  be  analyzed  at  supersonic  speeds. 

2.2  Types  of  Data  Computed 

2.2.1  Aerodynamics  -  The  program  computes  the  following  aerodynamic  parame¬ 
ters  as  a  function  of  angle  of  attack  for  each  configuration: 

Normal  Force  Coefficient 
Lift  Coefficient 

C  Pitching  Moment  Coefficient 

m 

X  Center  of  Pressure  In  calibers  from  the  moment  reference  center 

cp 

C.  Axial  Force  Coefficient 

A 

Cjj  Drag  Coefficient 

Cy  Side  Force  Coefficient 

C  Yawing  Moment  Coefficient 

n 

C^  Rolling  Moment  Coefficient 

C,.  Normal  force  coefficient  derivative  with  angle  of  attack 

Na 

C^  Pitching  moment  coefficient  derivative  with  angle  of  attack 

C„  Side  force  coefficient  derivative  with  sideslip  angle 

•H 

C  Yawing  moment  coefficient  derivative  with  sideslip  angle 

C.  Rolling  moment  coefficient  derivative  with  sideslip  angle 


5 


The  derivative  output  can  be  in  degrees  or  radians.  Partial  output  results, 
which  detail  the  components  used  in  the  calculations,  are  also  optionally 
available. 

The  program  has  the  capability  to  perform  a  static  trim  of  the  configura¬ 
tion,  using  any  finset  for  control  with  fixed  incidence  on  the  other  sets. 
The  two  types  of  aerodynamic  output  available  from  the  trim  option  are  as 
follows: 

0  Untrimmed  data  -  Each  of  the  aerodynamic  force  and  moment  coeffi¬ 
cients  are  printed  in  a  matrix,  which  is  a  function  of  angle  of 

attack  and  panel  deflection  angle.  This  output  is  optional. 

0  Trimmed  data  -  The  trimmed  aerodynamic  coefficients,  and  trim  deflec¬ 
tion  angle,  are  output  as  a  function  of  angle  of  attack. 

2.2.2  Geometry  -  All  components  of  the  configuration  have  their  physical  pro¬ 
perties  calculated  and  output  for  reference  if  requested.  All  data  is  sup¬ 

plied  in  the  user  selected  system  of  units. 

2.2.3  Other  -  The  reference  area  and  reference  length  are  user  defined.  The 
user  may  optionally  select  to  print  the  calculated  body  or  fin  pressure  coef¬ 
ficient  distributions  at  supersonic  speeds.  Outputs  of  the  partial  aerody¬ 
namic  results  and  a  summary  of  method  extrapolations  are  also  optionally 
available. 

2.3  Operational  Considerations 

The  code  has  been  written  to  conform  to  the  coding  standards  for  the  Amer¬ 
ican  National  Standards  Institute  (ANSI)  FORTRAN  IV  or  FORTRAN  66.  The  code 
has  also  been  designed  to  be  compatible  with  ANSI  Standard  FORTRAN  X3. 9-1978, 
often  referred  to  as  FORTRAN  V  or  FORTRAN  77. 
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There  are  only  four  exceptions  to  the  ANSI  standards  used  in  the  computer 
code: 

1.  Transfer  on  end-of-file.  The  FORTRAN  IV  statement  IF(EOF(UNIT) .NE.O) 
is  used  in  the  standard  COC  compatible  code.  The  FORTRAN  V  statement 
EN0=labe1  in  the  FORTRAN  READ  statements  is  incorporated  in  the  code,  but 
is  inactive  in  the  COC  code. 

2.  NAMELIST.  The  use  of  namelist  for  input  and  output  (I/O)  is  used. 
Although  this  appears  to  be  a  violation  of  FORTRAN  IV,  it  is  really  not 
since  a  namelist  emulator  has  been  written  for  the  Missile  Datcom  code 
using  FORTRAN  IV. 

3.  Mathematical  functions.  A  few  mathematical  functions  are  not  con¬ 
sidered  "standard,"  such  as  the  trigonometric  tangent.  Standard  FORTRAN 
equivalents  for  these  functions  are  available  on  request. 

4.  PROGRAM  card.  The  code  was  developed  on  the  Control  Data  Corporation 
CYBER  computers.  This  system  requires  that  the  first  card  of  the  main 
routine  be  a  PROGRAM  statement.  An  IBM  or  VAX  compatible  version  of  the 
code  is  also  maintained  which  has  a  different  format  for  the  program  card. 

Procedures  for  making  the  code  operational  on  other  than  Control  Data 
Corporation  computer  equipment  are  given  in  Volume  2  to  this  manual.  Only 
minor  coding  changes  are  required  and  can  be  accomplished  by  anyone  having  an 
understanding  of  the  FORTRAN  language.  The  code  is  currently  operational  on 
CYBER.  VAX.  PRIME,  IBM.  UNIVAC  and  HARRIS  computer  systems. 
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3.0  INPUT  DEFINITION 


Inputs  to  the  program  are  grouped  by  "case.“  A  "case"  consists  of  a  set 
of  input  cards  which  define  the  flight  conditions  and  geometry  to  be  run. 
Provisions  are  made  to  allow  multiple  cases  to  be  run.  The  successive  cases 
can  either  Incorporate  the  data  of  the  previous  case  (using  the  input  card 
SAVE)  or  be  a  completely  new  configuration  design.  The  SAVE  feature  permits 
the  user  to  define  a  body  and  wing  (or  canard)  configuration  in  the  first  case 
and  vary  the  tail  design  for  subsequent  cases. 

The  scheme  used  to  input  data  to  the  computer  program  is  a  mixture  of 

namelist  and  control  cards.  (See  Section  4.5  on  alternative  input  schemes.) 

This  combination  permits  the  following: 

0  Inputs  are  column  independent  and  can  be  input  in  any  order. 

0  All  numeric  Inputs  are  related  to  mnemonic  (variable)  names. 

0  Program  input  "flags"  are  greatly  reduced.  Required  "flags"  are 

identified  by  a  unique  alphabetic  name  which  corresponds  to  the  op¬ 
tion  selected. 

The  program  includes  an  error  checking  routine  which  scans  all  inputs  and 
identifies  all  errors.  This  process  is  a  single-pass  error  checking  routine; 
all  errors  are  identified  in  a  single  "run."  In  addition,  the  program  checks 
for  necessary  valid  Inputs,  such  as  a  non-zero  Reynolds  number.  In  some 
cases,  the  code  will  take  corrective  action.  The  type  of  corrective  action 
taken  is  summarized  later  in  this  section. 

Flexibility  has  been  maintained  for  all  user  inputs  and  outputs.  The 
following  summarize  the  program  generality  available: 

0  The  units  system  can  be  feet,  inches,  meters  or  centimeters.  The 

default  is  feet. 
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Derivatives  can  be  expressed  In  degree  or  radian  measure.  Degree 
measure  is  the  default. 


0  The  body  geometry  can  be  defined  either  by  shape  type  or  by  surface 
coordinates. 

0  The  airfoil  can  be  user  defined.  NACA  or  supersonic  shaped  sections. 
The  NACA  sections  are  defined  using  the  NACA  designation.  A  hexagon- 
ally  shaped  supersonic  section  Is  the  default. 

0  The  configuration  can  be  run  at  a  fixed  sideslip  angle  and  varying 

body  angle  of  attack,  or  a  fixed  aerodynamic  roll  angle  and  varying 
total  angle  of  attack. 

0  The  flight  conditions  can  be  user  defined,  or  set  using  a  Standard 

Atmosphere  model.  The  capability  to  define  wind  tunnel  test  condi¬ 
tions  as  the  flight  conditions  is  also  available.  The  default  flight 
condition  is  zero  altitude. 

3.1  Namelist  Inputs 

The  required  program  inputs  use  FORTRAN  namelist.  Missile  Oatcom  is  simi¬ 
lar  to  other  codes  which  use  the  namelist  Input  technique,  but  differ  as 
follows: 

0  Namelist  Inputs  are  column  Independent,  and  can  begin  In  any  column 
Including  the  first.  If  a  namelist  Is  continued  to  a  second  card, 
the  continued  card  must  leave  column  1  blank.  Also,  the  card  before 
the  continued  card  must  end  with  a  comma.  The  last  usable  column  Is 
number  79  if  column  1  Is  used,  and  column  80  If  column  1  Is  blank. 

0  The  same  namelist  can  be  Input  multiple  times  for  the  same  Input 
case.  The  total  number  of  namelists  read.  Including  repeat  occurr¬ 
ences  of  the  same  namelist  name,  must  not  exceed  300. 
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The  three  namelist  Inputs 


JREFQ  SREF=1.,$ 

$REFQ  LREF=2.,$ 

$REFQ  R0UGH=0.001.$ 

are  equivalent  to 

$REFQ  SREF=1 ..LREF-2.,R0U6H=0.001,$ 

The  last  occurence  of  a  namelist  variable  In  a  case  Is  the  value  of  that  vari¬ 
able  used  for  the  calculations. 

The  three  namelist  Inputs 

$REFQ  SREF=1.,$ 

$FLTC0N  NMACH=2.,MACH=1.0,  2.0,$ 

$REFQ  SREF=2..$ 

are  equivalent  to 

$REFQ  SREF=2.,$ 

$FLTC0N  NMACH=2.,  MACH=1.0,  2.0,$ 

0  The  namelists  can  be  Input  In  any  order. 

0  Only  those  namelists  required  to  execute  the  case  need  be  entered. 

0  Certain  hollerlth  constants  are  permitted.  They  are  summarized  In 
Table  3.  Note  that  any  variable  can  be  Initialized  by  using  the  con¬ 
stant  UNUSED;  for  example,  LREF=UNUSED  sets  the  reference  length  to 
Its  Initialized  value. 
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All  Missile  Oatcom  namelist  inputs  are  either  real  numbers  or  logical  con¬ 
stants.  Integer  constants  will  produce  a  nonfatal  error  message  from  the 
error  checking  routine  and  should  be  avoided. 

The  namelist  names  have  been  selected  to  be  mnemonically  related  to  their 
physical  meaning.  The  eight  namelists  available  are  as  follows: 


Namelist 

Inputs 

$FLTC0N 

Flight  Conditions  (Angles  of  attack,  Mach 

numbers,  etc.) 

)REFQ 

Reference  quantities  (Reference  area, 

length,  etc.) 

$AXIB0O 

Axisymmetric  body  definition 

$ELLB0D 

Elliptical  body  definition 

IFINSETn 

Fin  descriptions  by  fin  set  (n  is  the  fin 
set  number;  1 ,  2,  3  or  4) 

IDEFLCT 

Panel  incidence  (deflection)  values 

ITRIM 

Trimming  information 

$INLET 

Inlet  geometry 

$EXPR 

Experimental  data 

$ARBB0O 

Arbitrary  body  geometry 

Each  component  of  the 

configuration  requires  a  separate  namelist  input. 

Hence,  an  input  case  of  a 

body-wing-tall  configuration  requires  at  least  one 

of  each  of  the  following  namelist  inputs,  since  not  all  variables  have  default 
values  assigned: 

$FLTC0N 

to  define  the  flight  conditions 

SAX I 600,  ARBB0O 

or  SELLB0D 

to  define  the  body 

SFINSETl 

to  define  the  most  forward  fin  set 

SFINSET2 

to  define  the  first  following  fin  set 

SFINSET3 

to  define  the  second  following  fin  set 

SFINSET4 

to  define  the  third  following  fin  set 
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The  following  namelists  are  optional  since  defaults  exist  for  all  inputs: 


SREFQ 

to  define  the  reference  quantities 

SDEFLCT 

to  define  the  panel  incidence 

angles) 

STRIM 

to  define  a  trim  case 

SINLET 

to  define  inlet  geometry 

$EXPR 

to  define  experimental  input  data 

(deflection 


Defaults  for  all  namelists  should  be  checked  to  verify  the  configuration  being 
modeled  does  not  include  an  unexpected  characteristic  introduced  by  a  default. 

The  following  sections  describe  each  of  the  namelist  inputs.  Each  section 
is  accompanied  by  a  figure  which  summarizes  the  input  variables,  their  defini¬ 
tions,  and  units.  Since  the  system  of  units  can  be  optionally  selected,  the 
column  "Units"  specifies  the  generic  system  of  units  as  follows: 

1  Units  of  length;  feet,  inches,  centimeters 

or  meters 

F  Units  of  force;  pounds  or  Newtons 

deg  Units  of  degrees;  if  angular,  in  angular 

degrees;  if  temperature,  either  degrees  Ran- 
kine  or  degrees  Kelvin 

sec  Units  of  time  in  seconds 

2 

Exponents  are  added  to  modify  the  above.  For  example,  I  means  units  of 

length  squared,  or  area.  Combinations  of  the  above  are  also  used  to  specify 

2 

other  units.  For  example,  F/i  means  force  divided  by  area,  which  is  a 
pressure. 

Since  it  is  difficult  to  discern  the  difference  between  the  number  zero 
and  the  alphabetic  letter  "0,"  this  manual  will  always  express  the  alphabetic 
letter  as  "0"  In  general,  the  number  zero  and  the  letter  "0"  are  not 
interchangable  unless  so  stated. 
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The  program  ascertains  the  configuration  being  modeled  by  the  presence  of 
each  component  namelist,  even  If  no  data  Is  entered.  The  following  rules  for 
namelist  Input  apply: 

0  Do  not  Include  a  namelist  unless  It  Is  required.  Once  read,  the 
presence  of  a  namelist  (and,  hence,  a  configuration  component)  can 
only  be  removed  using  the  DELETE  control  card  In  a  subsequent  case. 
Simply  setting  all  variables  to  their  Initialized  values  will  not 
remove  the  configuration  component. 

0  Do  not  Include  a  variable  within  a  namelist  unless  It  Is  required. 
Program  actions  are  often  determined  from  the  number  and  types  of 
Input  provided. 

o  Do  not  over-spec  If  y  the  geometry.  User  Inputs  will  take  precedence 
over  program  calculations.  Inputs  that  define  a  shape  that  is  phy¬ 
sically  Impossible  will  be  used  as  specified.  The  program  does  not 
"fix-up"  inconsistent  or  contradictory  Inputs. 

3.1.1  NAMELIST  FLTC0N  -  Flight  Conditions 

This  namelist  defines  the  flight  conditions  to  be  run  for  the  case.  The 
program  Is  limited  to  no  more  than  20  angles  of  attack  and  20  Mach  numbers  per 
case  at  a  fixed  sideslip  angle,  aerodynamic  roll  angle,  altitude,  and  panel 
deflection  angle.  Therefore,  a  "case"  Is  defined  as  a  fixed  geometry  with 
variable  Mach  number  and  angles  of  attack. 

The  inputs  are  given  in  Figure  1.  There  are  two  ways  in  which  the  aero¬ 
dynamic  pitch  and  yaw  angles  can  be  defined: 

0  Input  ALPHA  and  BETA.  If  BETA  is  input  and  PHI  Is  not.  It  Is  assumed 
that  the  body  axis  angles  of  attack  (a)  and  sideslip  angles  (  )  are 
defined. 
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0  Input  ALPHA  and  PHI.  If  PHI  is  input  and  non-zero,  it  is  assumed 
that  ALPHA  is  the  total  angle  of  attack  (a)  and  PHI  is  the  aero¬ 
dynamic  roll  angle  . 

0  Input  ALPHA,  BETA  and  PHI.  The  value  for  BETA  is  ignored  if  PHI  is 
non-zero. 

As  a  minimum  the  following  variables  must  be  defined: 

NALPHA  number  of  angles  of  attack  to  run  (NALPHA>2) 

ALPHA  angle  of  attack  schedule  (matching  NALPHA) 

NNACH  number  of  Mach  numbers  or  speeds  (NMACH>1} 

MACH  or  VINE  Mach  number  or  speed  ^^chedule  (matching  NMACH) 

The  REN,  TINF  and  PINF  data  must  correspond  to  the  MACH  or  VINE  inputs.  The 
ALPHA  and  MACH  dependent  data  can  be  input  in  any  order;  the  code  will  sort 
the  data  into  ascending  order. 

Reynolds  number  is  always  required.  Three  types  of  inputs  are  permitted 
to  satisfy  the  Reynolds  number  requirement: 

(1)  Specify  Reynolds  number  per  unit  length  using  REN 

(2)  Specify  the  altitude  using  ALT,  and  the  speed  using  MACH  or  VINE 
(Reynolds  number  is  computed  using  the  Standard  Atmosphere  model) 

(3)  Specify  pressure  and  temperature  using  PINF  and  TINF,  and  the 
speed  using  MACH  or  VINE  (typical  of  data  available  from  a  wind 
tunnel  test) 

User  supplied  data  will  take  precedence  over  program  calculations.  Hence, 
the  user  can  override  any  default  or  Standard  Atmosphere  calculation.  The 
default  condition  is  sea-level  altitude  (ALT=0.)  if  the  wrong  combination  of 
inputs  are  provided  for  the  Reynolds  number  cannot  be  calculated. 
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3.1.2  NAMELIST  REFQ  -  Reference  Quantities 

Inputs  for  this  namelist  are  optional  and  are  defined  in  Figure  2.  A  ve¬ 
hicle  scale  factor  (SCALE)  permits  the  user  to  Input  a  geometry  that  is  scaled 
to  the  size  desired.  This  scale  factor  is  used  as  a  multiplier  to  the  user 
defined  geometry  inputs;  it  is  not  applied  to  the  user  input  reference  quanti¬ 
ties  (SKFF,  LREF,  LATREF).  If  no  reference  quantities  are  input,  they  are 
computed  based  upon  the  scaled  geometry.  XCG  is  input  relative  to  the  origin 
(X=0)  and  is  scaled  using  SCALE. 

In  lieu  of  specifying  the  surface  roughness  height  R0UGH,  the  surface 
Roughness  Height  Rating  (RHR)  can  be  specified.  The  RHR  represents  the  arith¬ 
metic  average  roughness  height  variation  in  millionths  of  an  inch.  Typical 
values  of  R0OUGH  and  RHR  are  given  in  Table  4. 

3.1.3  Namelist  AXIB0D  -  Axisymmetric  Body  Geometry 

An  axisymmetric  body  is  defined  using  this  namelist.  The  body  can  be 
specified  in  one  of  two  ways: 

OPTION  1 ;  The  geometry  is  divided  into  nose,  centerbody,  and  aft  body 
sections.  The  shape,  overall  length,  and  base  diameter  for  each  section 
are  specified.  Note  that  not  all  three  body  sections  need  to  exist  on  a 
configuration;  for  example,  a  nose-cylinder  configuration  does  not  require 
definition  of  an  aft  body. 

OPTION  2:  The  longitudinal  stations  and  corresponding  body  radii  are  de¬ 
fined,  from  nose  to  tail. 

It  is  highly  recommended  that  Option  1  be  used  when  possible.  The  program 
automatically  calculates  the  body  contour  based  upon  the  segment  shapes  using 
geometry  generators.  Hence,  more  accurate  calculations  are  possible.  Even 
when  Option  2  is  used,  appropriate  Option  1  inputs  should  be  included.  This 
identifies  where  the  code  should  insert  break  points  in  the  contour. 
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The  namelist  input  variables  are  given  in  Figure  3  and  a  sketch  of  the 
geometric  inputs  are  given  in  Figure  4.  The  program  uses  the  input  value  for 
NX  to  determine  which  option  is  being  used.  If  NX  is  not  input  then  Option  1 
inputs  are  assumed.  If  both  shapes  and  body  coordinates  (Options  1  and  2)  are 
used,  the  body  coordinate  information  will  take  precedence.  NX  can  be  set  to 
its  initialized  value  (to  simulate  the  variable  as  not  input)  by  specifying 
"NX=UNUSE0''. 

Most  of  the  subsonic  and  transonic  methods  are  dependent  upon  nose  and 
afterbody  size  and  shape.  Hence,  Option  2  inputs  should  also  specify  the  var¬ 
iables  LN0SE,  DN0SE,  LCENTR,  DCENTR,  LAFT,  and  DAFT  to  further  define  the 
geometry.  If  these  parameters  are  not  input,  they  are  selected  as  follows: 

LN0SE  -  Length  of  the  body  segment  to  where  the  radius  first  reaches 
a  maximum 

ON0SE  -  The  diameter  at  the  first  radius  maximum 

LCENTR  -  Length  of  the  body  segment  where  the  radius  is  constant 

OCENTR  -  Diameter  of  the  constant  radius  segment 

LAFT  -  The  remaining  body  length 

DAFT  -  Diameter  at  the  base 

DEXIT  -  Not  defined  (implies  that  base  drag  is  not  to  be  included  in 

the  axial  force  calculations) 

If  DEXIT  is  not  input,  or  set  to  UNUSED,  the  base  drag  computed  for  the 
body  geometry  will  not  be  included  In  the  final  computed  axial  force  calcula¬ 
tions.  To  Include  a  "full"  base  drag  Increment,  a  zero  exit  diameter  must  be 
specified  (DEXIT=0.). 

If  body  coordinates  are  input  using  the  variables  NX,  X,  R,  and  OISC0N, 
and  the  nose  is  spherically  blunted,  the  geometry  must  be  additionally  defined 
using  the  following: 

0  BN0SE  must  be  specified  (even  if  zero) 
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0  TRJNC  must  be  set  to  .FALSE. 

0  The  first  five  (5)  points  in  the  X  and  R  arrays  must  lie  on  the 

spherical  nose  cap  [i.e.,  X(l),  X(2),  X(3),  X(4),  X(5),  R(l),  R(2), 
R(3),  R(4),  and  R(5)  are  spherical  cap  coordinates] 

The  lollowing  summarizes  the  Input  generality  available: 

0  X(l)  does  not  have  to  be  0.0;  an  arbitrary  origin  can  be  selected. 

0  Five  shapes  can  be  specified  by  name: 

C0NICAL  (CONE)  -  cone  or  cone  frustrum  (default  for  boattails  and 
flares) 

0GIVE  -  tangent  ogive  (default  for  noses)  POWER  -  power  law* 

P0WER  -  power  law* 

HAACK  -  L-V  Haack  (length-volume  constrained)* 

KARMAN  -  von  Karman  (L-0  Haack;  length-diameter  constrained)* 

0  If  the  nose  is  truncated  (TRUNC=.TRUE. )  it  is  assumed  that  the  nose 
is  open  (i.e.,  the  forward  flat  face  produced  will  not  be  a  drag  pro¬ 
ducing  element) 

0  If  DAFT<0CENTR  the  afterbody  is  a  boattail.** 

0  If  0AFT>0CENTR  the  afterbody  Is  a  flare.** 

0  If  LAFT  is  not  input,  aft  body  (boattail  or  flare)  does  not  exist. 

*  applies  to  noses  only 

**  DAFT  must  not  be  equal  to  DCENTR 
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3.1.4  NAMELIST  ELLB0D  -  Define  Elliptically-Shaped  Bodies 

Elliptically-shaped  cross  section  bodies  are  defined  using  this  namelist. 
The  inputs  are  similar  to  those  for  the  axisymmetric  body  geometry  (AXIB0O), 
and  are  shown  in  Figure  5a-5b.  The  types  of  shapes  available,  and  the  limita¬ 
tions,  are  the  same  as  those  given  for  axisymmetric  bodies.  Please  read  Sec¬ 
tion  3.1.3  for  limitations. 

Note  that  the  body  cross  section  ellipticity  can  vary  along  the  body  long¬ 
itudinal  axis.  Sections  which  are  taller-than-wide  and  wider-than-tall  can  be 
mixed  to  produce  "shaped"  designs.  The  shape  of  the  sections  is  controlled  by 
the  variables  EN0SE,  ECENTR,  and  EAFT  or  ELLIP,  H  and  W. 

3.1.5  NAMELIST  FINSETn  -  Define  Fin  Set  n 

Figure  6a  describes  the  variables  needed  to  be  input  for  fin  set  planform 
geometry  descriptions.  Optional  fin  cross-section  inputs  are  described  in 
Figure  6b.  Special  user  specified  fin  cross-sections  can  be  input  using  the 
variables  in  Figure  6c.  The  user  may  specify  up  to  four  non-overlapping  fin 
sets.  The  variable  "n"  in  the  namelist  specifies  the  fin  set  number.  Fin 
sets  must  be  numbered  sequentially  from  the  front  to  the  back  of  the  mission 
beginning  with  fin  set  one.  An  input  error  will  occur  if  "n"  is  zero  or 
omitted.  The  code  allows  for  between  1  and  8  geometrically  identical  panels 
to  be  input  per  fin  set.  The  panels  may  be  arbitrarily  rolled  about  the  body 
and  can  be  given  dihedral. 

Four  types  of  airfoil  sections  are  permitted — hexagonal  (HEX),  circular 
arc  (ARC),  NACA  airfoils  (NACA)  and  user  defined  (USER).  Only  one  type  of 
airfoil  section  can  be  specified  per  fin  set,  and  this  type  is  used  for  all 
chord  wise  cross  sections  from  root  to  tip.  Diamond-shaped  sections  are  con¬ 
sidered  a  special  case  of  the  HEX  type;  hence,  hexagonal  and  diamond  sections 
can  coexist  on  the  same  panel.  The  airfoil  proportions  can  be  varied  from 
span  station  to  span  station. 
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The  user  selects  "break  points"  on  the  panel  (Figure  7).  A  "break  point" 
specifies  a  change  in  leading  or  trailing  edge  sweep  angle.  Also  a  break 
point  may  specify  a  change  in  airfoil  section,  but  the  section  must  be  of  the 
same  type  (i.e.,  a  change  in  section  type  cannot  go  from  a  NACA  to  an  ARC) 
only  the  proportions  can  change.  The  location  of  each  "break  point"  is  defin¬ 
ed  by  specifying  its  semi-span  station  (SSPAN)  from  the  vehicle  centerline  and 
position  of  the  chord  leading  edge  (XLE)  from  the  first  body  station.  The 
"break  point"  chord  leading  edge  array  (XLE)  can  be  defined  by  simply  specify¬ 
ing  the  root  chord  leading  edge  [XLE(l)]  and  the  sweep  angles  of  each  succes¬ 
sive  panel  segment  if  the  semi-span  stations  are  input.  Note  that  only  those 
variables  that  uniquely  define  the  fin  need  to  be  entered.  Redundant  inputs 
can  lead  to  numerical  inconsistencies  and  subsequent  computational  errors. 

The  panel  sweep  angle  (SWEEP)  can  be  specified  at  any  span  station  for 
each  segment  of  the  panels.  If  STA=0.,  the  sweep  angle  input  is  measured  at 
the  segment  leading  edge;  if  STA=1.,  the  sweep  angle  input  is  measured  at  the 
segment  trailing  edge.  Note  that  some  aerodynamic  methods  are  very  sensitive 
to  panel  sweep  angle.  For  small  span  fins,  small  errors  in  the  planform  in¬ 
puts  can  create  large  sweep  angle  calculation  errors.  It  is  recommended  that 
exact  sweep  angles  be  specified  wherever  possible;  for  example,  if  the  panel 
trailing  edge  is  unswept,  specifying  SWEEP=0.  and  STA=1 .  will  minimize  calcu¬ 
lation  error.  Then  the  leading  edge  sweep  will  be  computed  by  the  code  inter¬ 
nally  using  the  SSPAN  and  CHORD  inputs. 

Since  all  panels  are  assumed  to  be  planar  (i.e.,  no  tip  dihedral),  all 
inputs  must  be  "true  view."  Once  the  planform  of  a  single  panel  is  defined, 
all  fins  of  the  set  are  assumed  to  be  identical.  The  number  of  panels  present 
is  defined  using  the  variable  NPANEL.  Each  panel  may  be  rolled  to  an  arbi¬ 
trary  position  around  the  body  using  the  variable  PHIF.  PHIF  is  measured 
clockwise  from  top  vertical  center  (looking  forward  from  behind  the  missile) 
as  shown  in  Figure  8.  Each  panel  may  also  contain  a  constant  dihedral.  A 
panel  has  zero  dihedral  when  it  is  aligned  along  a  radial  ray  from  the  center- 
line  (see  Figure  10).  The  variable  used  to  specify  dihedral  is  GAM.  GAM  is 
positive  if  the  panel  tip  chord  is  rotated  clockwise  (see  Figure  10). 
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Different  aerodynamics  will  be  computed  depending  upon  whether  the 
FLTC0N  namelist  variable  PHI,  or  the  FINSETn  namelist  variable  PHIF,  is  used 
to  roll  the  geometry.  Figure  8  depicts  the  usage  of  the  roll  options.  The 
variable  "PHI"  means  that  the  body  axes  system  is  to  be  rolled  with  the  mis¬ 
sile  body,  whereas  PHIF  keeps  the  aerodynamics  in  a  non-rolled  body  axis,  but 
rather  locates  the  fin  positions  around  the  body.  PHIF  must  be  input  for  each 
panel,  while  PHI  rolls  the  whole  configuration. 

It  is  the  user's  responsibility  to  assure  that  the  fins  are  (1)  on  the 
body  surface,  and  (2)  do  not  lie  Internal  to  the  body  mold  line.  The  program 
does  not  check  for  these  peculiarities.  If  SSPAN(1)=0  is  input,  the  program 
will  assume  that  the  panel  semi-span  data  relative  to  Its  root  chord  are  sup¬ 
plied.  The  code  will  automatically  Interpolate  the  body  geometry  to  place  the 
panel  on  the  body  surface  with  the  root  chord  parallel  to  the  body  center- 
line.  See  Section  3.4  for  modeling  fins  on  body  segments  of  varying  radii. 

When  defining  more  than  one  fin  set,  the  fin  sets  must  never  have  their 
planforms  overlap  one  another.  There  must  be  sufficient  space  between  the 
forward  fin  trailing  edge  and  aft  fin  leading  edge  to  avoid  violating  the  as¬ 
sumptions  made  by  the  aerodynamic  computations.  It  is  assumed  by  the  aero¬ 
dynamic  model  that  the  vortices  are  fully  rolled  up  when  they  pass  the  control 
points  of  the  next  downstream  set  of  fins.  In  reality  the  vortex  sheet  does 
not  fully  roll  up  until  it  is  at  least  four  semispans  downstream.  If  two  fin 
sets  are  closer  than  this  the  results  may  be  in  error  since  the  use  of  a  vor¬ 
tex  filament  model  may  introduce  too  much  vorticity.  The  closer  the  spacing 
the  larger  the  error  may  be.  No  algorithm  error  will  result  from  too  close  a 
fin  set  spacing. 

Panels  with  cut-out  portions  can  be  modeled  by  using  one  of  the  ten  avail¬ 
able  fin  segments  as  a  transition  segment.  This  is  accomplished  by  giving  the 
segment  a  small  span,  such  as  0.0001,  and  specifying  the  segment  root  and  tip 
chords  to  transition  into  the  cut-out  portion  of  the  fin. 
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3.1.6  NAHELIST  DEFLCT  -  Panel  Deflection  Angles 

This  namelist  permits  the  user  to  fix  the  incidence  angle  for  each  panel 
in  each  fin  set.  The  variables  are  given  in  Figure  9.  Note  that  the  panel 
numbering  scheme  is  assumed  to  be  that  shown  in  Figure  10.  The  array  element 
of  each  deflection  array  corresponds  to  the  panel  number. 

The  scheme  for  specifying  deflection  angles  is  unique,  yet  concise.  The 
scheme  used  is  based  upon  the  body  axis  rolling  moment: 

"In  Missile  Datcom  a  positive  panel  deflection  is  one  which  will  produce  a 

negative  (counterclockwise  when  viewed  from  the  rear)  roll  moment  incre¬ 
ment  at  zero  angle  of  attack  and  sideslip." 

3.1.7  NAMELIST  TRIM  -  Trim  Aerodynamics 

This  namelist  instructs  the  program  to  statically  trim  the  vehicle  longi¬ 
tudinally  (C^=0)-  The  inputs  are  given  in  Figure  11.  Note  that  only  one 
fin  set  can  be  used  for  trimming.  The  user  only  specifies  the  range  of  de¬ 
flection  angles  desired  using  DELHIN  and  OELMAX;  the  code  will  try  to  trim  the 
vehicle  for  each  angle  of  attack  specified  using  the  allowable  fin  deflec¬ 
tions.  This  option  will  not  trim  the  vehicle  at  a  specific  angle  of  attack  if 
the  deflection  required  is  outside  the  range  set  by  the  values  of  OELMIN  and 
DELMAX. 

The  deflection  sign  convention  used  is  that  described  in  Section  3.1.6; 
hence,  OELMIN  and  DELMAX  are  input  as  if  deflecting  the  panel  to  the  maximum 
will  produce  a  negative  rolling  moment  from  the  panels  resulting  normal  force 
Increment.  The  magnitude  of  DELMIN  must  always  be  less  than  the  magnitude  of 
DELMAX. 

A  logical  variable,  ASYM,  has  been  included  to  permit  reverse  panel  de¬ 
flections.  For  example,  deflecting  all  panels  in  one  sense  results  in  a  roll¬ 
ing  moment  and  no  pitching  moment.  The  ASYM  flag  will  permit  analysis  of  an 
elevator  (or  pitch  deflection)  effect,  by  deflecting  panels  on  one  side  of 


the  vehicle  only,  with  opposite  panels  mirroring  those  deflections.  Since  a 
maximum  of  eight  panels  are  allowed  In  a  f Inset,  only  four  panels  of  the  fin- 
set  can  be  deflected  in  the  reverse  direction  using  the  ASYM  flag.  Both 
trimmed  and  untriramed  results  are  available  for  output. 

3.1.8  NAMELIST  INLET  -  Axisymmetric  and  2-D1mens1onal  Inlet  Geometry 

The  namelist  Inputs  shown  In  Figure  12  will  model  both  2-  dimensional  (2- 
0)  and  axisymmetric  Inlet  Installations  of  the  type  depicted  in  Figure  13, 
The  variables  INTYPE  and  POSO  are  used  to  Identify  the  type  of  Inlet  and  the 
attachment  scheme.  If  the  diverter  attaches  to  the  top  of  a  2-D  Inlet,  the 
Inlet  Is  a  top  attached  2-0  Inlet.  If  the  diverter  attaches  to  the  side  of 
the  2-0  Inlet,  the  Inlet  Is  a  side  attached  2-0  Inlet.  For  an  axisymmetric 
inlet,  this  distinction  Is  unnecessary  because  of  the  symmetry. 

Figure  14  Illustrates  the  procedure  Missile  Oatcom  uses  to  obtain  the  pro¬ 
per  Inlet  orientation  based  upon  the  Input  data.  To  provide  capability  for 
more  general  Inlet  geometries  to  be  Incorporated  in  the  future,  the  Inlet  in¬ 
stallation  Is  considered  to  consist  of  two  major  components:  an  Inlet  and  a 
diverter  which  attaches  the  Inlet  to  the  missile  body.  To  setup  the  namelist 
Inputs,  the  Inlet  and  diverter  are  assumed  to  be  mounted  In  the  horizontal 
position  on  the  missile  body  as  shown  at  the  top  of  Figure  14.  Missile  Datcom 
will  combine  the  diverter  and  Inlet  inputs  Into  a  fully  Integrated  Inlet  In¬ 
stallation  at  this  position,  and  then  rotate  1t  to  the  specified  angular 
orientation,  PHI. 

The  variable  NIN  defines  the  number  of  Inlets  being  modeled.  The  program 
can  handle  1,  2,  3,  and  4  Inlets  In  combination.  Figure  15  illustrates  the 
type  PHI  values  for  single,  twin,  triple  and  quadruple  inlet  configurations. 
For  four  Inlet  combinations,  a  positive  angle  should  be  input  for  PHI(l)  and  a 
negative  angle  for  PHI(2)  to  obtain  cruciform  inlets.  For  Inlet  geometries 
with  inlets  in  the  vertical  plane,  an  angle  of  90^  should  be  used;  the  pro¬ 
gram  will  then  position  one  inlet  on  the  bottom  vertical  and  other  on  the  top 
vertical . 
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The  namelist  inputs  of  inlet  coordinates  and  dimensions  require  the  defi¬ 
nition  of  an  inlet  coordinate  system  (Xj,  Y^,  Z^)  and  its  location  rela¬ 
tive  to  the  missile  coordinate  system  (X  ,  Y  ,  Z  ).  Figure  16  shows  the 
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coordinate  systems  for  the  three  inlet  installations  assuming  they  are  ini¬ 
tially  in  the  horizontal  position.  For  all  inlets,  the  inlet  coordinate  sys¬ 
tem  axial  location  is  located  relative  to  X=0  with  input  X  ,  (Although  the 

BU 

origin  shown  at  the  inlet  leading  edge,  it  can  be  at  any  arbitrary  X-sta- 

tion).  For  the  side  attached  2-0  Inlet  the  origin  location  relative  to  the 

missile  horizontal  centerline  is  located  by  coordinate  Z__.  The  coordinates 
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Yj  arid  Zj  relative  to  the  missile  are  computed  from  the  input  geometry. 

2-D  Inlet  -  The  nomenclature  for  the  2-0  inlet  is  provided  in  Figure  17.  The 
inlet  is  considered  to  be  composed  of  four  major  sections:  forebody,  midbody, 
boattail  and  base.  The  user  identifies  these  sections  through  inputs  of  se¬ 
lected  coordinates  and  dimensions  at  locations  1  through  10.  The  forebody  is 
the  most  complicated  section  and  includes  the  top  and  bottom  lip  leading  edge 
which  form  the  top  and  bottom  of  the  inlet  opening.  The  external  surface  on 
the  top  and  outside  edge  can  be  chamfered.  (The  side  external  transfer  can 
end  at  or  before  the  inlet  maximum  depth  location).  The  maximum  depth  of  the 
inlet  occurs  on  the  forebody.  Internal  to  the  forebody  at  the  top  of  the  in¬ 
let  opening  is  the  internal  compression  ramp. 

The  midbody  is  a  rectangular  tube  connecting  the  forebody  and  boattail. 
The  boattail  begins  with  a  rectangular  cross  section  and  transitions  into  a 
semicircular  base  section. 

Figure  18  shows  a  side  attached  2-0  Inlet,  the  19  inputs  required  to  de¬ 
fine  its  geometry,  and  the  set  of  assumptions  used  to  construct  the  inlet  ge¬ 
ometry.  For  example,  the  forebody  and  midbody  sides  are  assumed  parallel. 
Therefore,  the  width  of  the  midbody  top  at  the  start  of  the  boattail  is  B^. 
Because  the  midbody  top  and  bottom  are  assumed  parallel,  the  coordinates  at 

location  3  and  8  are  Z_=Z-,  X  -X_  and  Z-=Z_.  The  leading  edge  of 
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the  internal  compression  ramp  begins  at  location  1,  the  top  lip  leading  edge, 
and  is  assumed  parallel  to  the  Xj  axis  to  location  9,  i.e.,  Zg*Z.j .  The 
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compression  ramp  is  assumed  to  terminate  at  the  axial  location  of  the  bottom 
lip  and,  therefore  X,-=X,.. 

I U  3 

In  the  program,  the  diverter  is  Integrated  with  the  inlet  and  missile 
body.  Figure  19  defines  the  side  attached  diverter  nomenclature,  inputs  and 
assumptions.  The  diverter  leading  edge  width,  TIND^ ,  sets  the  spacing  be¬ 
tween  the  body  and  inlet  side  and  XINO^  sets  the  diverter  leading  edge  axial 

location  in  the  inlet  coordinate  system.  The  vertical  location  is  assumed  to 

be  at  the  half  height  of  the  inlet  midbody.  The  only  other  input  required  is 
the  axial  position,  XIND^.  at  which  the  diverter  height  becomes  equal  to  the 
inlet  midbody  height.  The  diverter  height  is  assumed  constant  from  XIND^  to 
the  start  of  the  inlet  boattail  at  X^.  Because  of  this  assumption,  XIND^ 
must  not  exceed  X^.  The  diverter  then  fairs  into  the  boattail  and  has  a 

height  at  the  inlet  base  equal  to  the  base  diameter.  Other  diverter  parameter 
such  as  width  are  computed  based  upon  these  assumptions. 

Top  Attached  Inlet  -  Figure  20  shows  the  18  inputs  required  to  define  the  ge¬ 
ometry  of  the  top  attached  2-0  inlet.  This  is  an  inlet  with  the  diverter 
attached  to  its  top  surface.  The  base  semicircle  is  oriented  with  the  flat 
side  on  top  to  provide  a  flat  top  surface  to  Integrate  with  the  diverter. 

The  diverter  construction  is  somewhat  different  than  for  the  side-mounted 

inlet.  Figure  21  defines  the  nomenclature.  Inputs  and  assumptions.  The  three 

diverter  Inputs  serve  the  same  function  as  for  the  side  mounted  2-D  inlet. 

The  other  dimensions  are  computing  assuming  their  inlet  midbody  is  parallel 

with  the  missile  centerline.  The  diverter  maximum  height  is  set  equal  to  the 

midbody  height,  B.,  and  its  base  height  Is  set  to  2’*B..  The  input  value  of 

2  4 

XIND^  must  not  be  less  than  the  Inlet  coordinate  X^  and  XIND^  must  not 
be  greater  than  X.. 

Axisvmmetric  Inlet  -  Figure  22  describes  the  nomenclature  for  the  axisymmetric 
inlet.  Five  major  inlet  sections  arc  defined:  centerbody,  lip,  midbody, 
boattail,  and  base.  The  centerbody  and  lip  are  modeled  by  inputs  at  3  axial 
stations  to  provide  an  approximate  description  of  their  actual  contours. 
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Figure  23  shows  the  15  inputs  required  to  define  the  geometry  and  the 
assumptions  made  internal  to  the  program  to  construct  the  remaining  geometry. 

through  B,  and  represent  radii  at  the  location  indicated.  B.  is 
the  offset  of  the  base  from  the  midbody  X.  through  X_  represent  the  iden- 
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tified  axial  stations.  Both  the  centerbody  and  cowl  lip  are  modeled  as  coni¬ 
cal  frustums.  The  midbody  of  the  inlet  is  assumed  cylindrical.  Therefore, 

the  radicS  as  station  7  is  equal  to  B  . 

b 

The  diverter  is  integrated  with  the  inlet  and  missile  body  as  shown  in 

Figure  24.  The  diverter  inputs  include  the  leading  edge  width,  XIN0,j  which 

sets  the  spacing  between  the  body  and  inlet  side,  and  XIN0,| ,  which  sets  the 

diverter  leading  edge  axial  location  in  the  inlet  coordinate  system  (XIN0,j 

must  be  between  X,  and  X,  so  that  the  diverter  starts  on  the  midbody). 

The  only  other  input  required  is  the  axial  position,  XINO^,  at  which  the 

diverter  height  becomes  equal  to  the  inlet  midbody  diameter  which  is  2B  . 

b 

the  geometry  for  the  remainder  of  the  diverter  is  constructed  assuming  a  con¬ 
stant  diverter  height  from  station  XINO^  until  the  start  of  the  inlet  boat- 
tail,  At  the  inlet  base  the  diverter  height  is  assumed  to  equal  2Bg.  Other 
diverter  parameters  and  widths  are  computed  based  upon  these  assumptions. 

3.1,9  NAMELIST  ARBB0D  -  Arbitrary  Body  Geometry 

Namelist  ARBB0D  allows  the  user  more  flexibility  in  creating  nonconven- 
tional  missile  body  shapes.  Body  shapes  input  should  not  be  too  flat  and 
should  have  fineness  ratios  greater  than  five  and  cross-sectional  aspect 
ratios  less  than  five. 


Geometric  reference  values  (e.g. 


etc)  and  aerody- 
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namic  conditions  (a,  M<»,  Re,  etc)  are  input  using  the  REFQ  and  FLTCON 
NAMELISTS.  These  NAMELISTs  have  already  been  described.  The  user  is  directed 
to  their  description  for  details  of  their  input  variables. 


The  geometry  of  a  configuration  is  developed  using  the  variables  described 
in  Figure  25.  Although  it  is  not  possible  to  describe  every  arbitrary  shape 
with  the  present  method,  it  does  expand  the  class  of  shapes  that  can  be  ana¬ 
lyzed  by  NAMELISTs  AXIB0O  and  ELLB0O. 
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The  planform  and  forebody  base  cross-section  are  the  most  important  shape 
details  needed  for  this  method.  The  planform  must  be  trapezoidal  or  triangu¬ 
lar,  to  simplify  the  input.  Consequently,  the  nose  must  converge  either  to  a 
point  or  to  a  flat  spatular  tip.  No  nose  bluntness  is  allowed  for  either  the 
pointed  or  spatular  nose  geometries. 

The  base  cross-section  of  the  nose  is  defined  by  connecting  one  or  more 
linear  segments  and/or  elliptic  arc  segments.  The  start  and  end  points  of  the 
segment  define  the  extent  of  that  segment.  Subsequent  segments  require  only 
the  end  point,  since  the  segment  begins  at  the  point  where  the  previous  seg¬ 
ment  ends. 

Mirror  symmetry  about  the  vertical  plane  reduces  the  chance  of  erroneous 
input  data.  Therefore,  only  that  portion  of  the  base  with  zero  or  positive  y 
values  may  be  input.  That  portion  which  falls  to  the  left  of  the  vertical 
plane  will  be  generated  by  the  program.  At  the  present  time,  only  eight  seg¬ 
ments  may  be  used  to  describe  the  half  cross-section.  The  first  segment  must 
start  at  the  x-z  plane,  proceed  counterclockwise,  and  end  at  the  x-z  plane. 

The  centerbody  must  have  the  same  cross-section  as  the  base  of  the  fore¬ 
body.  The  only  input  required  to  describe  the  centerbody  therefore  is  the 
length.  No  afterbody  is  permitted  at  the  present  time. 

To  illustrate  the  inputs  required  to  describe  the  arbitrary  cross  section 
bodies,  four  examples  are  presented  next.  The  first  example  is  very  simple. 
The  second  illustrates  a  pointed,  arbitrary  cross  section  nose.  The  third 
details  the  inputs  for  a  spatular  nose  and  the  fourth  describes  another  form 
of  a  pointed,  arbitrary  cross  section  nose. 

The  forebody  in  Figure  26  has  a  pointed  tip,  and  a  circular  base  cross 
section.  The  forebody  length  is  1.111  ft  and  the  centerbody  length  is  2.222 
ft,  The  tip  width  is  zero  and  is  placed  at  the  origin.  The  nose  is  a  right 
circular  cone,  so  the  center  of  the  base  lies  on  the  x-axis;  the  radius  is 
0.1667  ft. 
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The  segment  may  be  described  in  two  ways.  First,  the  angle  of  the  radius 
at  the  starting  point  is  -90  degrees,  relative  to  the  positive  x-axis;  the  end 
point  is  +90  degrees  from  the  x-axis.  Likewise,  the  y-z  coordinate  pair  for 
the  starting  point  is  (0.0000,  -  .1667),  and  (0.0000,  +  .1667)  for  the  end 
point. 

Only  one  segment  is  used  to  describe  the  base,  and  the  base  type  for  that 
segment  is  3.  The  tip  is  a  point,  so  the  tip  type  must  be  zero.  The  calcula¬ 
tion  of  C  requires  the  integration  of  along  the  longitudinal 
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axis  segment  by  segment.  The  code  divides  the  planform  into  equal  length  seg¬ 
ments  for  which  C.,  is  calculated. 

The  forebody  in  Figure  27  also  has  a  pointed  tip,  but  the  cross  section  is 
a  rounded  equilateral  triangular.  Except  for  the  base  radius,  which  is  not 
applicable  here,  the  reference  lengths  are  the  same  as  in  Figure  1. 

Four  segments,  two  linear  and  two  elliptic  arcs,  are  used  to  describe  the 
basic  half  cross-section.  The  first  segment  is  linear,  beginning  at  (0.0000, 
-  .1700)  and  extending  to  (0.1200,  -  .1700).  The  second  segment  is  a  circular 
arc;  the  origin  is  located  at  (0.1200,  -.1033).  The  arc  extends  from  -90 
degrees  to  +30  degrees.  The  third  segment  is  a  line,  extending  from  the  end 
of  the  second  segment  to  (0.0578,  0.1367).  The  fourth  segment  is  a  circular 
arc  with  a  radius  of  0.0667.  The  end  point  may  be  specified  either  as  +90 
degrees,  or  as  (0.0000,  0.1700).  The  base  types  are  2,  3,  2  and  3  respective¬ 
ly;  the  tip  is  a  point  so  all  four  tip  types  must  be  zero. 

The  third  example.  Figure  28,  is  a  wedge  with  a  trapezoid  planform  and 
base  cross-section.  The  forebody  length  is  0.6460,  and  because  there  is  no 
centerbody,  that  length  is  zero.  The  nose  tip  is  spatular;  the  total  tip 
width  is  0.1667  ft.  The  center  of  the  tip  has  been  placed  at  the  origin,  so 
the  distance  to  the  tip  is  zero. 

The  base  half  cross-section  may  be  described  using  three  line  segments. 
The  first  starts  at  (0.0000,  0.0000)  and  extends  to  (0.2088,  0.0000).  The 
second  segments  continues  from  the  end  of  the  first  segment  to  (0.0833, 
0.0565).  The  third  segment  ends  on  the  z-axis  at  (0.0000,  0.0565). 
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Since  all  three  segments  are  linear,  all  three  base  types  are  2.  The 
first  segment  lies  In  the  x-y  plane,  and  narrows  to  a  finite  width  at  the  tip, 
so  the  tip  type  is  1.  The  second  segment  slopes  up,  and  to  the  left,  and  nar¬ 
rows  to  a  point,  so  the  2nd  tip  type  is  0.  the  final  segment  narrows  to  a 
finite  width  at  the  tip,  so  the  third  tip  type  is  1. 

Figure  29  shows  a  forebody  with  a  pointed  tip  and  a  half  ellipsoid  base. 
The  forebody  is  1.111  ft  and  the  centerbody  is  2.222  ft  long. 

Two  segments  describe  the  base  half  cross-section;  an  elliptsoid  and  a 
linear  segment.  Both  extend  to  a  point  at  the  nose  the  tip  types  are  both 
zero.  The  linear  segment  starts  at  (0.000,  0.000)  and  goes  to  (.2886, 
0.000).  The  ellipsoid  segment  automatically  begins  at  the  end  of  the  last 
segment,  therefore  no  OSTART  value  must  be  specified.  The  ellipsoid  segment 
is  centered  around  (0.000,  0.000)  and  has  major  and  minor  axes  of  lengths 
.2886  and  .1924.  The  ellipsoid  is  traced  from  0  through  90  degrees,  therefore 
0FINAL=90,O. 

3.1.10  NAMELIST  EXPR  -  Experimental  Data  Substitution 

This  namelist  is  used  to  substitute  experimental  data  for  the  theoretical 
data  generated  by  the  program.  The  variables  to  be  input  ae  shown  in  Figure 
30.  Use  of  namelist  EXP  does  not  stop  the  program  from  calculating  theoreti¬ 
cal  data,  but  rather  the  experimental  data  is  used  in  configuration  synthesis, 
and  it  is  the  experimental  data  that  is  used  for  the  component  aerodynamics 
for  which  it  is  input. 

Experimental  data  may  be  substituted  for  any  configuration  component  or 
partial  configuration.  Experimental  data  is  input  at  a  specific  Mach  number, 
when  using  namelist  EXPR  the  case  must  be  run  at  the  Mach  number  for  which  you 
are  substituting  experimental  data.  However,  the  experimental  data  being  in¬ 
put  may  have  different  reference  quantities  and  a  different  center  of  gravity 
location  than  the  case  being  run. 
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Experimental  data  input  for  a  fin  alone  is  input  as  panel  data  -  not  as 

total  fin  set  data.  The  user  should  note  that  experimental  data  for  fin  alone 

C  is  not  used  in  the  configuration  synthesis  process.  Instead  fin  alone 
ma 

C.,  (the  experimental  value  if  input)  i  used  to  determine  the  fin  contri- 

Na 

bution  to  C  during  configuration  synthesis.  If  body  alone  experimental 

rfkk 

data  and  body-fin  experimental  data  are  input  for  the  same  case  the  body  data 

is  ignor-^d  in  configuration  synthesis.  If  experimental  C  data  is  input 

ma 

for  a  body  +  1  fin  set  for  a  multi-fin  set  configuration,  the  calculated  con¬ 
tributions  to  C  of  the  other  fin  sets  are  added  to  the  experimental  data, 
ma 

Since  the  experimental  namelist  forms  the  basis  for  configuration  incre¬ 
menting.  the  lateral  directional  coefficients  are  included  to  allow  for  side¬ 
slip  cases.  These  coefficients  are  input  the  same  as  the  longitudinal  coeffi¬ 
cients.  However,  if  the  lateral  directional  coefficients  are  input,  the 
lateral  directional  beta  derivatives  will  not  be  computed  our  output. 

The  following  rules  apply  to  the  use  of  namelist  EXPR. 

1.  It  is  assumed  that  the  coefficients  in  EXPR  are  for  the  same  sideslip 
and/or  aerodynamic  roll  as  the  case  being  run. 

2.  Separate  namelist  EXPR  must  be  specified  for  each  Mach  number. 

3.  Each  namelist  EXPR  must  end  with  a  $END  card. 

4.  Separate  namelist  EXPR  must  be  specified  for  each  partial  configura¬ 
tion  for  which  experimental  data  is  to  be  input,  (i.e.,  body,  body  + 
1  fin  set,  etc)  since 

5.  Separate  namelist  EXPR  must  be  specified  for  each  reference  quantity 
change. 

Example: 

The  user  has  experimental  data  available  for  a  body  +  2  fin  set  con¬ 
figurations  and  is  Interested  in  the  effects  of  adding  a  booster  con¬ 
taining  a  third  fin  set.  he  would  then  use  namelist  EXPR  to  input 
the  experimental  data.  When  the  configuration  is  synthesized,  it 
would  use  the  experimental  data  for  body  +  2  fin  sets  and  theoretical 
data  for  fin  set  three. 
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3.2  Control  Card  Inputs 

Control  cards  are  one  line  commands  which  select  program  options.  Al¬ 
though  they  are  not  required  inputs,  they  permit  user  control  over  program 
execution  and  the  types  of  output  desired.  Control  cards  enable  the  following: 

0  Printing  internal  data  array  results  for  diagnostic  purposes  (DUMP) 

0  Outputting  intermediate  calculations  (PART,  BUILD,  PRESSURES,  PRINT 

AER0.  PRINT  EXTRAP,  PRINT  6E0M,  PL0T,  NAMELIST,  WRITE,  F0RMAT) 

0  Selecting  the  system  of  units  to  be  used  (DIM,  DERIV) 

0  Defining  multiple  cases,  permitting  the  reuse  of  previously  input 
namelist  data  or  deleting  namelists  of  a  prior  case  (SAVE,  DELETE, 
NEXT  CASE) 

0  Adding  case  titles  or  comments  to  the  input  file  and  output  pages 
(*,  CASEID) 

0  Limits  the  calculations  to  longitudinal  aerodynamics  (NO  LAT) 

3.2.1  Control  Card  -  General  Remarks 

A  total  of  42  different  control  cards  are  available.  There  is  no  limit  to 
the  number  of  control  cards  that  can  be  present  in  a  case.  If  two  or  more 
control  cards  contradict  each  other,  the  last  control  card  input  will  take 
precedence.  All  control  cards  must  be  input  as  shown,  including  any  blanks. 
Control  cards  can  start  in  any  column  but  they  cannot  be  continued  to  a  second 
card.  Misspelled  cards  are  ignored.  Control  cards  can  be  located  anywhere 
within  a  case. 

Once  input,  the  following  control  cards  remain  in  effect  for  all  subse¬ 
quent  cases: 
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DIM  FT 

DIM  IN 

DIM  CM 

DIM  M 

FORMAT 

HYPER 

INCREMENT 

N060 

NO  LAT 

PLOT 

SOSE 

WRITE 

The 

following  control  cards  are 

effective 

only  for  the  case 

in  which  they 

appear: 

BUILD 

CASEID 

DAMP 

DELETE 

DUMP  CASE 

DUMP  NAME 

NAMELIST 

PART 

PRESSURES 

PRINT  AERO 

PRINT  EXTRAP 

PRINT  GEOM 

SAVE 

SPIN 

TRIM 

The 

control  cards 

can  be  changed 

from  case 

to  case: 

DERIV  DEG 

DERIV  RAD 

NACA 

The  only  control  card  that  can  be  optionally  saved,  from  case-to-case,  is 
the  NACA  card. 

Available  control  cards  are  summarized  as  follows: 

BUILD 


This  control  card  Instructs  the  program  to  print  the  results  of  a  configu¬ 
ration  build-up.  All  configurations  which  can  be  built  from  the  components 
defined  will  be  synthesized  and  output,  including  Isolated  data  (e.g.,  body 
alone,  fin  alone,  etc.).  Component  build-up  data  is  not  provided  if  the  TRIM 
option  is  selected. 

CASEID 


A  user  supplied  title  to  be  printed  on  each  output  page  is  specified.  Up 
to  72  characters  can  be  specified  (card  columns  8  to  80). 
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DAMP 


When  DAMP  control  card  is  input  longitudinal  dynamic  derivatives  are  com¬ 
puted  and  the  results  output  for  the  configuration.  Dynamic  derivatives  are 
configuration  components  or  partial  configurations  may  be  output  using  the 
PART  or  BUILD  control  cards  respectively. 

DELETE  namel.nameZ 

This  control  card  instructs  the  program  to  ignore  a  previous  case  namelist 
input  that  was  retained  using  the  SAVE  control  card.  All  previously  saved 
namelists  with  the  names  specified  will  be  purged  from  the  input  file.  Any 
new  inputs  of  the  same  namelist  will  be  retained.  At  least  one  name  (namel) 
must  be  specified. 

DIM  IN.  DIM  FT.  DIM  CM  or  DIM  H 

This  control  card  sets  the  system  of  units  for  the  user  inputs  and  program 
outputs.  The  four  options  are  inches  (DIM  IN),  feet  (DIM  FT),  centimeters 
(DIM  CM),  and  meters  (DIM  M) .  The  default  system  of  units  is  feet.  Once  the 
system  of  units  has  been  set,  it  remains  set  for  all  subsequent  cases  of  the 
"run". 

DERIV  DEG  or  DERIV  RAD 

All  output  derivatives  are  set  to  either  degree  (DERIV  DEG)  or  radian 
(DERIV  RAD)  measure*.  The  derivative  units  can  be  changed  more  than  once  dur¬ 
ing  the  run  by  inputting  multiple  DERIV  cards. 

DUMP  CASE 


Internal  data  blocks,  used  in  the  computation  of  the  case,  are  written  on 
Tape  6.  This  control  card  automatically  selects  partial  output  (PART). 

*  Default  is  degree. 
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DUMP  namel  .nanie2 


This  permits  the  user  to  write  selected  internal  data  blocks  on  Tape  6. 
At  least  one  name  (namel)  must  be  specified.  Dump  names  are  identified  in 
Volume  2  -  Program  Implementation  Guide.  The  arrays  will  be  dumped  in  units 
of  feet,  pounds,  degrees  or  degrees  Rankine. 

F0RMAT  (format) 

This  control  card  is  used  in  conjunction  with  the  WRITE  control  card.  It 
specifies  the  format  of  the  data  to  be  printed  to  tape  unit  3.  The  format  is 
input  starting  with  a  left  parenthesis,  the  format  and  a  right  parenthesis. 
This  is  exactly  the  same  as  a  FORTRAN  FORMAT  statement.  Because  of  the  code 
structure  alphanumeric  data  must  not  be  printed.  For  example 

F0RMAT  ((8(2X.F10.4))  is  legal 
F0RMAT  (2HX=,F10.4)  is  illegal 

The  default  format  is  8F10.4,  and  will  be  used  if  the  F0RMAT  control  card  is 
not  present.  Multiple  formats  can  be  used.  The  last  F0RMAT  read  will  be 
used  for  all  successive  WRITE  statements  until  another  F0RMAT  is  encounter¬ 
ed.  Hence,  the  F0RMAT  must  precede  the  applicable  WRITE. 

HYPER 


This  control  card  causes  the  program  to  select  the  Newtonian  flow  method 
for  bodies  at  any  Mach  number  above  1.4.  HYPER  should  normally  be  selected  at 
Mach  numbers  greater  than  6. 

INCRMT 


This  card  is  used  to  set  the  configuration  Incrementing  flag.  Configura¬ 
tion  incrementing  uses  the  first  case  of  a  run  to  determine  correction  factors 
for  the  longitudinal  and  lateral  aerodynamic  coefficients.  These  correction 
factors  are  computed  by  comparing  theoretical  and  experimental  values  for  each 
coefficient  for  which  data  is  input.  The  experimental  values  are  input  using 
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el1st  EXPR.  During  subsequent  cases  of  the  run,  the  correction  factors  are 
applied  to  coefficients  for  which  experimental  data  was  input  in  the  first 
case.  This  provides  the  user  with  a  method  to  evaluate  changes  in  a  configu¬ 
ration. 

The  INCRMT  card  must  be  input  in  the  first  case  of  a  run.  The  first  case 
must  be  run  at  the  same  Mach  number  as  the  experimental  data  which  is  input. 
Once  the  increment  flag  is  set  it  can  not  be  deleted  during  that  run. 

The  following  restrictions  apply: 

1.  All  cases  of  a  run  must  have  the  same  number  of  fin  sets. 

2.  All  cases  of  a  run  must  have  the  same  sideslip  or  aerodynamic  roll 
angle  as  the  first  case  (BETA  or  PHI  as  specified  in  namelist  FLTCON). 

3.  The  first  case  must  be  run  at  exactly  the  same  angles  of  attack  as 
the  experimental  data  being  input. 

4.  All  cases  must  be  run  within  the  same  Mach  regime  (subsonic,  tran¬ 
sonic,  or  supersonic)  as  the  experimental  data. 

5.  Experimental  data  can  only  be  input  in  the  first  case  and  only  for 
the  complete  configuration.  No  additional  data  can  be  substituted. 

6.  To  increment  C„  and  experimental  data  must  be  input  for 

Yrt  Na 

Cy  and  C^. 

Use  of  configuration  incrementing  may  or  may  not  increase  the  accuracy  of 
the  results.  The  following  guidelines  will  produce  better  results  when  using 
configuration  incrementing: 

1.  The  user  may  run  different  angles  of  attack  in  each  case.  However, 
no  angle  of  attack  should  exceed  the  upper  or  lower  limit  of  the 
angles  of  attack  for  which  experimental  data  was  input  in  the  first 
case. 

2.  Experimental  data  should  be  input  at  as  many  angles  of  attack  ?s 
possible. 
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3.  The  user  should  remember  that  the  effect  of  a  change  in  Mach  number 
from  case  to  case  is  not  corrected  by  inputting  experimental  data  at 
one  Mach  number  as  is  required. 

NACA 

This  card  defines  the  NACA  airfoil  section  designation  (or  supersonic  air¬ 
foil  definition).  Note  that  if  airfoil  coordinates  and  the  NACA  card  are 
specified  for  the  same  aerodynamic  surface,  the  airfoil  coordinate  specifica¬ 
tion  will  be  used.  Therefore,  if  coordinates  have  been  specified  in  a  previ¬ 
ous  case  and  the  SAVE  option  is  in  effect,  the  saved  namelist  must  be  deleted 
or  the  namelist  variable  SECTYP  must  be  changed  for  the  NACA  card  to  be  recog¬ 
nized  for  that  aerodynamic  surface.  The  airfoil  designated  with  this  card 
will  be  used  for  all  segments  and  panels  of  the  fin  set. 

The  form  of  this  control  card  and  the  required  parameters  are  as  follows: 


Card  Column(s) 

InDut(s) 

Purpose 

1  thru  4 

NACA 

The  unique  letters  NACA  designate 
that  an  airfoil  is  to  be  defined 

5 

Any  delimeter 

6 

1  or  2 

Fin  set  number  for  which  the  air¬ 
foil  designation  applies 

7 

Any  delimeter 

8 

1.4,5,6,S 

Type  of  NACA  airfoil  section; 
1-series  (1),  4-dig1t  (4),  5-digit 
6-ser1es  (6),  or  supersonic  (S) 

(5). 

9 

Any  delimeter 

10  thru  80 

Designation 

Input  designation  (see  Table  5); 
columns  are  free-field  (blanks  are 
nored) 

ig- 

Only  fifteen  (15)  characters  are  accepted  in  the  airfoil  designation.  The 
vocabulary  consists  of  the  following  characters: 

0123456789  A,-.- 
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Any  characters  input  that  are  not  in  the  vocabulary  list  will  be  interpreted 
as  the  number  zero  (0).  Table  5  details  the  restrictions  on  the  NACA  designa¬ 
tion. 

NAMELIST 


This  control  card  instructs  the  program  to  print  all  namelist  data.  This 
is  useful  when  multiple  inputs  of  the  same  variable  or  namelist  are  used. 

NEXT  CASE 


This  card  indicates  termination  of  the  case  input  data  and  instructs  the 
program  to  begin  case  execution.  It  is  required  for  multiple  case  "runs." 
This  card  must  be  the  last  card  input  for  the  case. 

N0  G0 

This  control  card  permits  the  program  to  cycle  through  all  of  the  input 
cases  without  computing  configuration  aerodynamics.  It  can  be  present  anywhere 
in  the  input  stream  and  only  needs  to  appear  once.  This  option  is  useful  for 
performing  error  checking  to  insure  all  cases  have  been  correctly  set  up. 

N0  LAT 

This  control  card  inhibits  the  calculation  of  the  lateral -directional 
derivatives  due  to  sideslip  angle.  Savings  in  computation  time  can  be  realiz¬ 
ed  by  using  this  option.  This  option  is  automatically  selected  when  using 
TRIM. 

PART 


This  control  card  permits  printing  of  partial  aerodynamic  output,  such  as 
a  summary  of  the  normal  force  and  axial  force  contributors.  Partial  output  of 
the  configuration  synthesis  methods  is  only  provided  if  the  TRIM  option  is  not 
selected.  Use  of  this  card  is  equivalent  to  inputing  all  PRINT  AERO  and  PRINT 
6E0M  control  cards  (except  *,  -,  +). 
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PL0T 


A  data  file  for  use  with  a  post-processing  plotting  program  is  provided. 
The  BUILD  control  card  must  also  be  present  in  oroer  to  generate  a  data  file 
with  component  buildup  data.  If  the  TRIM  option  is  selected,  only  the  trimmed 
and  untrimmed  total  configuration  results  are  provided. 

PRESSURES 


This  control  card  instructs  the  program  to  print  the  body  and  fin  alone 
pressure  coefficient  distributions  at  supersonic  speeds.  Only  pressure  data 
to  15  degrees  angle  of  attack  for  bodies  and  at  zero  angle  of  attack  for  fins 
are  printed. 

PRINT  AER0  name 


This  control  card  instructs  the  program  to  print  the  incremental  aerodyna¬ 
mics  for  "name,"  which  can  be  one  of  the  following: 


B0DY 

for  body  aerodynamics 

FINl 

for  FINSETl  aerodynamics 

FIN2 

for  FINSET2  aerodynamics 

FIN3 

for  FINSET3  aerodynamics 

FIN4 

for  FINSET4  aerodynamics 

SYNTHS 

for  configuration  synthesis  aerodynamics 

TRIM 

for  trim/untrimmed  aerodynamics 

BEND 

for  panel  bending  moments 

HINGE 

for  panel  hinge  moments 

INLET 

for  inlet  aerodynamics 

All  options  are  automatically  selected  when  the  control  card  PART  is  used. 
Details  of  the  output  obtained  with  these  options  are  presented  in  Section  4.2. 
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PRINT  EXTRAP 


This  control  card  enables  the  printing  of  method  extrapolation  messages 
produced  during  execution  of  the  case. Extrapolation  messages  are  not  normally 
provided. 

PRINT  Ge0M  name 

This  control  card  instructs  the  program  to  print  the  geometric  character¬ 
istics  of  the  configuration  component  "name,"  which  can  be  one  of  the  follow¬ 
ing: 


B0DY 

for  body  geometry 

FINl 

for 

FINSETl  geometry 

FIN2 

for 

FINSET2  geometry 

FIN3 

for 

FINSET3  geometry 

FIN4 

for 

FINSET4  geometry 

INLET 

for 

inlet  geometry 

All  options  are  automatically  selected  when  the  control  card  PART  is  used. 
SAVE 


The  SAVE  card  saves  namelist  inputs  from  one  case  to  the  following  case 
but  not  for  the  entire  run.  This  permits  the  user  to  build-up  or  change  a  com¬ 
plex  configuration,  case-to-case,  by  adding  new  namelist  cards  without  having 
to  re-input  namelist  cards  of  the  previous  case.  When  changing  a  namelist 
that  has  been  saved,  the  namelist  must  first  be  deleted  using  the  delete  con¬ 
trol  card. 


39 


The  only  control  card  that  can  be  optionally  saved,  case-to-case,  is  the 
NACA  card. 

S0SE 

The  presence  of  this  control  card  selects  the  Second-Order  Shock  Expansion 
Method  for  axl symmetric  bodies  at  supersonic  speeds.  S0SE  should  be  select¬ 
ed  if  any  Mach  number  is  higher  than  2.0. 

SPIN 

When  the  SPIN  control  card  Is  input,  spin  and  magnus  derivatives  are  com¬ 
puted  for  body  alone.  If  the  configuration  being  run  Is  a  body  +  fin  sets, 
the  spin  derivatives  are  still  computed  for  body  alone.  A  PART  or  BUILD  card 
must  be  input  for  body  alone  derivatives  to  be  printed  out. 

TRIM 

This  control  card  causes  the  program  to  perform  a  trim  calculation.  Com¬ 
ponent  buildup  data  cannot  be  dumped  If  TRIM  Is  selected.  The  use  of  this 
control  card  is  the  same  as  If  the  namelist  TRIM  was  Included  except  that  the 
defaults  for  namelist  TRIM  are  used. 

WRITE  name,  start,  end 

This  control  card  causes  the  common  block  "name"  to  be  printed  to  tape 
unit  3  using  the  most  recent  F0RMAT  control  card.  Locations  from  "start"  to 
"end"  are  dumped  (see  Volume  II  for  a  list  of  block  names).  Multiple  WRITE 
statements  may  be  Input,  and  there  Is  no  limit  to  the  number  which  may  be  pre¬ 
sent.  The  presence  of  a  WRITE  will  cause  the  block  "name"  to  be  printed  for 
all  cases  of  the  run.  The  output  will  be  In  units  of  feet,  pounds,  degrees, 
or  degrees  Rankine.  If  the  PL0T  option  Is  also  selected,  this  output  will 
be  "mixed"  with  the  PL0T  file  output  on  tape  unit  3. 

«r 

Any  card  with  an  asterisk  (*)  in  Column  1  will  be  Interpreted  as  a  comment 

card.  This  permits  detailed  documentation  of  case  inputs. 
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Figure  31  schematically  shows  how  Missile  Datcom  inputs  are  structured. 
This  example  illustrates  a  multiple  case  job  in  which  case  2  uses  part  of  the 
case  1  inputs.  This  is  accomplished  through  use  of  the  SAVE  control  card. 
Case  1  is  a  body-wing-tail  configuration;  partial  output,  component  buildup 
data,  and  a  plot  file  are  created.  Case  2  uses  the  body  and  tail  data  of  case 
1  (the  wing  is  deleted  using  DELETE),  specifies  panel  deflection  angles  and 
sets  the  data  required  to  trim. 

There  is  no  limit  to  the  number  of  cases  that  can  be  ''stacked"  In  a  single 
run,  provided  that  no  more  than  300  namelist  inputs  are  "saved"  between 
cases.  If  a  SAVE  control  card  is  not  present  in  a  case,  all  previous  case 
inputs  are  deleted. 

3.3.1  Configuration  Incrementing  Case  Setup 

A  "configuration  incrementing"  case  setup  is  shown  in  Figure  32.  This 
figure  shows  the  Inputs  for  a  three  case  setup  fin  parametric  analysis.  The 
first  case  is  the  calibration  case  with  the  remaining  cases  being  used  for  the 
parametric  analysis.  Therefore,  the  first  case  must  contain  both  the  INCRMT 
control  card  and  EXPR  namelist.  These  should  only  appear  in  the  first  case. 

3.4  Special  Usage  of  Input  Parameters 

It  is  possible  to  manipulate  the  input  geometry,  such  that  unique  configu¬ 
rations  can  be  modeled.  This  section  defines  those  capabilities. 

3.4.1  Locating  Panels  on  Varying  Body  Radii  Segments 

The  fin  panels  can  be  located  anywhere  on  the  geometry.  If  they  are  to  be 
positioned  on  a  varying  radii  segment,  select  the  root  chord  span  station 
(SSPAN(l)]  such  that  the  center  of  the  exposed  root  chord  is  on  the  surface 
mold  line.  Physically  this  places  part  of  the  panel  within  the  body  and  part 
offset  from  the  body. 
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If  SSPAN(1)  1s  input  precisely  as  zero,  the  code  will  assume  that  panel 
semi-span  stations  relative  to  the  root  chord  are  defined.  It  will  then  in¬ 
terpolate  the  body  geometry  at  the  root  chord  center  and  add  the  body  radius 
at  this  point  to  the  user  defined  values  in  the  SSPAN  array. 
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Table  1  Body  Addressable  Configurations 


M1 1-1 4065 


1 .  Nose  Shape 

Conical 
Sharp 
Blunted 
Tnrncated 
Tangent  Ogive 
Sharp 
Blunted 
Tmncated 
Other 

2.  Centerbody  Shape 

Cyfnder 

Bliptically  Variable 
Arbitrary  Cross  Section 

3.  Afterbody  Shape 


Conical 

Tangent  Ogive 
Other 
Flares 
Conical 
Ogive 
Other 
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Table  2  Subsonic/Transonic  Method  Limitations 


M1 1-14066 


Range  Permitted 

Action  Taken  If  Exceeded 

Nose  Shape 

CONE  or  OGIVE 

Uses  OGIVE 

Conical  Nose  Bluntness  (Ca) 

Sharp  Only 

Uses  Sharp  Method 

Nose  Bluntness  (0^,0^) 

Sharp  Only 

Uses  Sharp  Method 

Conical  Nose  Slope 

0  to  25  Degrees 

Uses  25  Degrees 

Boattail  Shape 

CONE  or  OGIVE 

Uses  CONE 

Conical  Boattail  Slope 

Oto  16  Degrees 

Extrapolates 

Ogive  Boatttail  Slope 

0  to  28  Degrees 

Extrapolates 

Flare  Shape 

CONE 

Uses  CONE 

Raid  Slope 

Oto  10  Degrees 

Extrapolates 

Airfoil  t/c 

Oto  12% 

Continues,  if  possibie 

Table  3  Namelist  Alphanumeric  Constants 


M11-14067 


NAME  LIST 

PERMITTED 

ALPHANUMERIC 

CONSTANTS 

CONVERTED 

VALUE 

(ALL) 

UNUSED 

1  .E-30  (Inttialized  Value) 

REFQ 

TURB 

0. 

NATURAL 

1. 

AXIB0D  or 

C0NICAL 

0. 

ELLB0D 

C0NE 

0. 

0QIVE 

1. 

P0WER 

2. 

HAACK 

3. 

KARMAN 

4. 

FINSETn 

HEX 

0. 

NACA 

1. 

ARC 

2. 

USER 

3. 

INLET 

AXI 

0. 

2D 

1. 
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Table  4  Equivalent  Sand  Roughness 


Ml  1-14068 


TYPE  OF  SURFACE 

EQUIVALENT  SAND 
ROUGHNESS 
k  (inches) 

RHR 

Aerodynamicalfy  smooth 

0 

0 

Polished  metal  or  wood 

0.02-0.08  X  10‘® 

6-26 

Natural  sheet  metal 

0.16  X  10*® 

53 

SiTxx)th  matte  paint,  carefully  applied 

0.25  X  10"® 

83 

Standard  camouflage  paint,  average  application 

0.40  X  10"® 

133 

Canrxiuflage  paint,  mass-production  spray 

1.20  X  10"® 

400 

Dip-galvanized  metal  surface 

6x10*® 

2000 

Natural  surface  of  cast  iron 

lOxlO'® 

3333 

Preferred  RHR  Values 


Application 

RHR 

Steel  structural  parts 

250 

Aluminum  an  titanium  stnjctural  parts 

125 

Close  tolerance  surfaces 

63 

Seals 

32 

VARIABLE  ARRAY 
NAME  DIMENSION 


NALPHA 

ALPHA 

BETA 

PHI 

NMACH 

MACH  ^ 

REN 

ALT 

VINF 

TINF 

PINF  > 


NAMELIST  FLTCON 


DEFINITION 


NUMBER  OF  ANGLES  OF  ATTACK  (AT  LEAST  2) 

ANGLE  OF  ATTACK  OR  TOTAL  ANGLE  OF  ATTACK 

SIDESLIP  ANGLE 

AERODYNAMIC  ROLL  ANGLE 

NUMBER  OF  MACH  NUMBERS  (AT  LEAST  1) 

MACH  NUMBERS 

REYNOLDS  NUMBER  PER  FOOT  (METER) 
GEOMETRIC  ALTITUDE 
FR^ESTREAM  SPEED 
FREESTREAM  STATIC  TEMPERATURE 
FREESTREAM  STATIC  PRESSURE 


Ml  1-14069 


UNITS  DEFAULT 


A  @ 
^/SEC@ 
DEG 
F/X2  0 


0  FOLLOWING  COMBINATIONS  SATISFY  THE  REYNOLDS  NUMBER  AND  MACH  NUMBER 
INPUT  REQUIREMENTS. 


USER  INPUT 

1. MACH,  REN 

2.  MACH,  ALT 

3.  VINF,  ALT 

4.  VINF,  TINF,  PINF 

5.  MACH,  TINF.  PINF 


(NONE) 

PINF,  TINF,  REN 
PINF,  TINF,  MACH,  REN 
MACH.  REN 
VINF.  REN 


©  INPUT  AS  1/FT  FOR  ENGLISH  UNITS  AND  1/M  FOR  METRIC  UNITS 
©  INPUT  AS  FT  FOR  ENGLISH  UNITS  AND  M  FOR  METRIC  UNITS 
0  INPUT  AS  FT/SEC  FOR  ENGLISH  UNITS  AND  M/SEC  FOR  METRIC  UNITS 
©  INPUT  AS  L»FT  2  for  ENGLISH  UNITS  AND  N/M  2  FOR  METRIC  UNITS 

Figure  1  Flight  Condition  Inputs 
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NAMELIST  REFQ 


M1 1-14070 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

SREF 

REFERENCE  AREA 

0 

© 

LREF 

— 

REFERENCE  LENGTH  (LONGITUDINAL) 

k 

LATREF 

- 

REFERENCE  LENGTH  (LATERAL-DIRECTIONAL) 

i® 

*LREF 

ROUGH 

RHR 

- 

SURFACE  ROUGHNESS  HEIGHT 

0. 

— 

ROUGHNESS  HEIGHT  RATING 

- 

0. 

XCG 

— 

LONGITUDINAL  POSITION  OF  C.G.  (+  AFT) 

k 

0. 

ZCG 

- 

VERTICAL  POSITION  OF  C.G.  (+  UP) 

k 

0. 

SCALE 

VEHICLE  SCALE  FACTOR  (MULTIPLIER  TO 
GEOMETRY) 

1. 

BLAYER 

BOUNDARY  LAYER  TYPE:  TURB  FOR  FULLY 
TURBULENT,  NATURAL  FOR  NATURAL  TRANSITION 

TURB 

©  EITHER  CAN  BE  USED. 


©  DEFAULT  IS  BODY  MAXIMUM  CROSS-SECTIONAL  AREA.  IF  NO  BODY  IS  INPUT, 
MAXIMUM  FIN  PANEL  AREA  IS  USED. 

©  DEFALUT  IS  BODY  MAXIMUM  DIAMETER.  IF  NO  BODY  IS  INPUT,  MAXIMUM 
FIN  PANEL  MEAN  GEOMETRIC  CHORD  IS  USED. 

©  INPUT  AS  INCHES  FOR  ENGLISH  UNITS  AND  CENTIMETERS  FOR  METRIC  UNITS. 


Figure  2  Reference  Quantity  Inputs 
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NAMELIST  AXIBOD 

OPTION  1 1NPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOORX0 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

0. 

TN0SE 

- 

NOSE  SHAPE  NAME: 

OGIVE 

CONICAL,  CONE,  OGIVE,  POWER, 

HAACK,  KARMAN 

P0WER 

- 

EXPONENT,  n,  FOR  POWER  SERIES  SHAPES, 

— 

0. 

(ri/R)  =  (Xi/A)n 

LN0SE 

- 

NOSE  LENGTH 

L 

DN0SE 

- 

NOSE  SECTION  BASE  DIAMETER 

A, 

1.0 

BN0SE 

- 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 

k 

0. 

NOSE  RADIUS 

TRUNC 

- 

TRUE.  IF  NOSE  IS  TRUNCATED 

— 

.FALSE 

LCENTR 

- 

CENTERBODY  LENGTH 

L 

0. 

DCENTR 

- 

CENTERBODY  BASE  DIAMETER 

k 

-  DNOSE 

TAFT  > 

AFTRBODY  SHAPE  NAME: 

— 

CONICAL 

1  ^ 

— 

CONICAL,  OGIVE 

LAFT  ^© 

- 

AFTERBODY  LENGTH 

k 

0. 

DAFT 

- 

AFTERBODY  BASE  DIAMETER 

k 

DExrr ) 

DIAMETER  OF  NOZZLE  EXIT 

k 

©  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  {i.e.  DAFT  ^  DCENTR) 


NAMELIST  AXIBOD 

OPTION  2  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOORX0 

LONnr,  UDINAL  COORDINATE  AT  NOSE  TIP 

k 

0. 

BN0SE 

— 

NOSt  BLUNTNESS  RADIUS  OR  TRUNCATED 

k 

0. 

- 

NOSE  RADIUS 

.FALSE 

TRUNC 

- 

TRUE.  IF  NOSE  IS  TRUNCATED 

— 

DExrr 

DIAMETER  OF  NOZZLE  EXIT 

k 

NX 

NUMBER  OF  INPUT  STATIONS 

X  © 

50 

LONGITUDINAL  COORDINATES 

k 

R 

50 

RADIUS  AT  EACH  X  STATION 

k 

DISC0N 

20 

INDICES  OF  X  STATIONS  WHERE  THE  SURFACE 

— 

SLOPE  IS  DISCONTINUOUS 

LNOSE  "10 

NOSE  LENGTH 

k 

LCENTR 

CENTERBODY  LENGTH 

k 

©  X  (NX)  MUST  BE  END  OF  BODY 


©  THESE  INPUTS  CAN  BE  USED  TO  IMPROVE  MODELLING  ACCURACY 


Figure  3  Axlsymmetrlc  Body  Geometry  Inputs 
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NAMELIST  ELLBOD 

OPTION  1 1NPUTS  M11  14072 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOORX0 

LONGITUDINAL  COORDINATE  NOSE  TIP 

0. 

TN0SE 

- 

NOSE  SHAPE  NAME: 

- 

OGIVE 

CONICAL  CONE  OGIVE  POWER 

HAACK  KARMAN 

POWER 

- 

EXPONENT,  n.  FOR  POWER  SERIES  SHAPES. 

- 

0. 

(ri/R)  =  (Xi/A)n 

LN0SE 

- 

NOSE  LENGTH 

k 

WN0SE 

- 

NOSE  SECTION  BASE  WIDTH 

k 

1. 

BN0SE 

— 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 

NOSE  RADIUS 

k 

0. 

TRUNC 

- 

-TRUE.  IF  NOSE  IS  TRUNCATED 

- 

•  FALSE 

EN0SE 

- 

ELLIPTICITY  AT  NOSE  BASE,  H/W 

- 

1. 

LCENTR 

- 

CENTERBODY  LENGTH 

k 

0. 

WCENTR 

CENTERBODY  BASE  WIDTH 

k 

=  WN0SE 

ECENTR 

— 

ELLIPTICITY  AT  CENTERBODY  BASE,  HW 

- 

1. 

TAFT  'I 

— 

AFTERBODY  SHAPE  NAME: 

- 

CONICAL 

— 

CONICAL  CONE  OGIVE 

LAFT 

AFTERBODY  LENGTH 

k 

0. 

WAFT 

AFTERBODY  BASE  WIDTH 

k 

EAFT  J 

ELLIPTICITY  AT  AFT  BODY  BASE,  YWJ 

- 

1. 

DEXIT 

DIAMETER  OF  NOZZLE  EXIT 

k 

Q  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  (i.e.  WAFT  ^  DCENTR) 


Figure  Sa  Elliptical  Body  Seometry  Inputs-Option  1 
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Ml  1-14073 


NAMELIST  ELLBOD 
OFnON  2  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOORX0 

- 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

B 

0. 

BN0SE 

- 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED  NOSE  RADIUS 

0. 

TRUNC 

- 

.  TRUE.  IF  NOSE  IS  TRUNCATED 

B 

.FALSE. 

NX 

- 

NUMBER  OF  INPUT  STATIONS 

■ 

X 

50 

LONGITUDINAL  COORDINATES 

w  O 

50 

BODY  HALF-WIDTH  AT  EACH  X  STATION 

B 

DISC0N 

20 

INDICES  OF  X  STATIONS  WHERE  THE  SURFACE  SLOPE 

IS  DISCONTINUOUS 

B 

H  O 

50 

BODY  HALF-HEIGHT  AT  EACH  X  STATION 

ELLIP  © 

50 

BODY  HEIGHT  TO  WIDTH  RATIO  AT  EACH  X  STATION 

B 

1. 

O  ONE  OF  THE  FOLLOWING  COMBINATIONS  IS  REQUIRED: 


W  AND  H.  W  AND  ELLIP,  OR  H  AND  ELLIP 


Figure  5b  Elliptical  Body  Geometry  Inputs-Optlon  2 


52 


Ml  1-1 4074 


NAMELIST  FINSETn 

NOMINAL  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

SECTYP 

— 

TYPE  OF  SECTION  TO  BE  DEFINED: 

HEX,  NACA,  ARC,  USER 

— 

HEX 

SSPAN  0 

10 

SEMI-SPAN 

k 

CHORD 

10 

PANEL  CHORD  LENGTH  AT  EACH  SSPAN 

k 

XLE 

10 

AXIAL  STATION  (X)  OF  CHORD  LEADINGL  EDGE 

AT  EACH  SSPAN 

k 

0.0 

SWEEP  0 

10 

SWEEPBACK  ANGLE  AT  EACH  SSPAN 

DEG 

0.0 

STA 

10 

CHORD  STATION  USED  IN  MEASURING  SWEEP  AT 

EACH  SSPAN  (0.  .  LEADING  EDGE,  1.  -  TRAILING  EDGE) 

— 

0.0 

LER0 

10 

PANEL  LEADING  EDGE  RADIUS  AT  EACH  SPAN  STATION 

k 

NPANEL 

— 

NUMBER  OF  PANELS  IN  SET  (1  -  8) 

— 

4 

PHIF0 

8 

ROLL  ANGLE  OF  EACH  FIN  MEASURED  CLOCKWISE 

FROM  TOP  VERTICAL  CENTER 

DEG 

— 

GAM 

8 

DIHEDRAL  OF  EACH  FIN  (POSITIVE  WHEN  PHIF  IS 
INCREASED,  SEE  FIG.  10) 

DEG 

0.0 

SKEW 

ANGLE  BETWEEN  THE  Y  AXIS  AND  THE  FIN  HINGE  LINE 
(POSITIVE  SWEPT  BACK) 

DEG 

0.0 

0  IF  SSPAN  (1)  -  0  INPUTS  ARE  RELATIVE  TO  ROOT  CHORD  NOT  BODY  CENTERLINE. 

0  IF  USING  SWEEP,  SPECIFY  ONLY  XLE  (1);  IF  USING  XLE  DO  NOT  SPECIFY  SWEEP. 

0  NOT  REQUIRED  FOR  NACA  AIRFOILS,  REQUIRED  FOR  USER  AIRFOILS 
0  IF  NOT  INPUT  THE  NUMBER  OF  PANELS  SPECIFED  ARE  EVENLY  SPACED  ABOUT  THE  BODY. 

Figure  6a  Fin  Geometry  Inputs-Moralnal 
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NAMELIST  FINSETn  Mii-i4075 


OPTIONAL  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

OEFINITiON 

UNITS 

DEFAULT 

ZUPPER 

10 

MAXIMUM  THICKNESS  OF  UPPER  SECTION  RATKDED 

TO  LOCAL  CHORD  LENGTH 

— 

0.025 

ZL0WER 

10 

MAXIMUM  THICKNESS  OF  LOWER  SECTION  RATKDED 
TO  LOCAL  CHORD  LENGTH 

— 

-ZUPPER 

LMAXU 

10 

FRACTION  CHORD  mOM  SECTION  LEADING  EDGE  TO 
MAXIMUM  THICKNESS  FOR  UPPER  SURFACE 

— 

0.5 

LMAXL 

10 

FRACTION  CHORD  FRCM  SECTION  LEADING  EDGE  TO 
MAXIMUM  THICKNESS  FOR  LOWER  SURFACE 

— 

-LMAXU 

LFLATU 

10 

FRACTION  CHORD  LENGTH  OF  CONSTANT  THICKNESS 
FOR  UPPER  SURFACE 

— 

0.0 

LFUTL 

10 

FRACTION  CHORD  LENGTH  OF  CONSTANTTHICKNESS 
FOR  LOWER  SURFACE 

— 

-LFUTU 

Figure  6b  Fin  Geometry  Inputs-Optlonal 
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NAMELIST  FlNSETn 

INPUTS  FOR  "USER"  SECTIONS  M11-14076 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XC0RD 

50 

PANEL  CHORD  STATION 

- 

MEAN  © 

50 

DISTANCE  BETWEEN  THE  MEAN  LINE  AND 
CHORD  LINE  AT  EACH  XCORD  FOR  AVERAGE 
AIRFOIL  SECTION 

THICK  © 

50 

THICKNESS  TO  LOCAL  CHORD  FRACTION  AT 
EACH  XCORD  FOR  AIRFOIL  SECTION 

- 

YUPPER  © 

50 

AIRFOIL  SECTION  UPPER  SURFACE 
COORDINATES,  FRACTION  CHORD,  AT  EACH 
XCORD 

— 

YLOWER  © 

50 

AVERAGE  AIRFOIL  SECTION  LOWER  SURFACE 
COORDINATES.  FRACTION  CHORD,  AT  EACH 
XCORD 

— 

NOTE:  ALL  VARIABLES  ARE  EXPRESSED  AS  UNITS  OF  CHORD 
O  EITHER  MEAN  AND  THICK  OR  YUPPER  AND  YLOWER  ARE  REQUIRED 
@  LEADING  EDGE  RADIUS  (VARIABLE  LER)  MUST  BE  DEFINED 


Figure  6c  Fin  Geometry  Input  User  Airfoils 
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view  Lo(rt(lno  Forward 
Figure  8  Roll  Attitude.  PHIF 


5i 


M1 1-1 4077 


NAMELIST  DEFLCT 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

DELTA1  0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  1  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

DELTA2  0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FM  SET  2  (SUBSCRIPT  IS  FW  NUMBER) 

DEG 

0. 

DELTAS  0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  3  (SUBSCRIPT  IS  FM  NUMBER) 

DEG 

0. 

DELTA4  0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FM  SET  4  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

XHINGE 

4 

AXIAL  STATION  (FROM  X  -  0)  OF  PANEL 

LINE 

k 

CR/2@ 

SKEW 

4 

SWEEPBACK  OF  HINGE  LINE 

DEG 

0. 

O  PANEL  #1  IS  TOP  (WHEN  <1>  -  0) 

^  DEFAULT  IS  AT  ONE-HALF  OF  THE  EXPOSED  ROOT  CHORD 
NOTE:  A  POSITIVE  DEFLECTION  ANGLE  PRODUCES  A  NEGATIVE  BODY 
AXIS  ROLLING  MOMENT  ATSRO  ANGLE  OF  ATTACK. 

Figure  9  Panel  Deflection  Inputs 


57 


z 


6W  (2)  •  0 


PANEL  2 

VIEW  LOOKING  FORWARD 


Figure  10  Fin  Numbering  and  Orientation 


NAMELIST  TRIM 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

SET 

— 

FIN  SET  TO  BE  USED  FOR  TRIM 

— 

1. 

PANL1  © 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

PANL2© 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

PANL3  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

PANL4  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

FALSE. 

PANL5  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

PANL6  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

PANL7  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

-FALSE. 

PANL8  0 

— 

.TRUE.  IF  PANEL  TO  BE  USED 

— 

.FALSE. 

DELMIN  0 

— 

MINIMUM  NEGATIVE  DEFLECTION 

DEG 

-25. 

DELMAX  0 

— 

MAXIMUM  POSITIVE  DEFLECTION 

DEG 

+20. 

ASYM 

4 

.TRUE.  IF  PANEL  WITH  SUBSCRIPT  IS  TO 

BE  DEFLECTED  OPPOSITE  TO  NORMAL  SIGN 
CONVENTION  (ASYMMETRIC  DEFLECTIONS) 

* 

.FALSE. 

©  DEFAULTS  APPLY  ONLY  IF  ALL  PANLX  DATA  ARE  NOT  INPUT  OR  .FALSE. 
©  BOTRDELMINANDDELMAXMUSTBESPECIFED 


Figure  11  Trim  Inputs 
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NAMELIST  INLET 


M1 1-14078 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

INTYPE 

— 

INLET  TYPE; 

2D  -  TWO  DIMENSIONAL  .  1 

AXI  -AXISYMMETRIC-0 

— 

0 

P0SD 

—  — 

FLAG  TO  IDENTIFY  THE  TYPE  OF  DIVERTER 

0  -  SIDE  MOUNTED  DIVERTER 

1  -  TOP  MOUNTED  DIVERTER 

—  — 

0. 

PHI(1)@ 

— 

ANGULAR  ORIENTATION  OF  INLETS  IN  SET  1 

DEG. 

0. 

PHI  (2)  @ 

— 

ANGULAR  ORIENTATION  OF  NLETS  IN  SET  2 

DEG. 

0. 

NIN 

— 

NUMBER  OF  INLETS  1. 2. 3.  OR  4 

— 

2. 

XB0OR 

XBO 

— 

LONGITUDINAL  DIST.  FROM  X  -  0  TO  INLET 
COORDINATE  SYSTEM  ORIGIN 

k 

— 

ZB0OR 

ZBO 

— 

DISTANCE  FROM  VEHICLE  CENTERLINE  TO 

INLET  COORDINATE  SYSTEM  ORIGIN 

k 

— 

© 

X 

10 

INLET  LONGITUDINAL  POSITIONS 

RELATIVE  TO  XBO 

k 

— 

B@ 

10 

INLET  RADIUS  OR  HALF  WIDTH 

k 

— 

Z@@ 

10 

INLET  HEIGHT  REUUVE  TO  ZBO 

k 

— 

XIND(1) 

—  ” 

DIVERTER  LEADING  EDGE  LOCATION 

RELATIVE  TO  INLET  COORDINATE 

SYSTEM  ORIGIN 

k 

—  — 

XIND  (2) 

— 

DIVERTER  MAXIMUM  HEIGHT  LOCATION 

RELATIVE  TO  INLET  COORDINATE 

SYSTEM  ORIGIN 

k 

— 

TIND 

— 

DIVERTER  LEADING  EDGE  WIDTH 

k 

— 

©  THE  VARIABLES  XI .  X2.  X3...  XI 0  CAN  BE  USED  INSTEAD  OF  ARRAY  ELEMENT  DESIGNATIONS. 
@  THE  VARIABLES  B1,  B2,  B3...  BIO  CAN  BE  USED  INSTEAD  OF  ARRAY  ELEMENT  DESIGNATIONS. 
©  THE  VARIABLES  Z1,  Z2.  Z3...  Z10  CAN  BE  USED  INSTEAD  OF  ARRAY  ELEMENT  DESIGNATIONS. 

§  MEASURED  POSITIVE  DOWNWARD  FROM  HORIZONTAL  POSITION. 

Z(1)  MUST  EQUAL  ZERO. 

Figure  12  Inlet  Geometry  Inputs 
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TOP  ATTACHED  2D  INLET 


SIDE  AHACHED  2D  INLET 


AXISYW1ETRIC  INLET 


Figure  13  Typical  Missile  Inlet  Installations 
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Z30  computed  by  code 


CONFIGURATION 

NIN 

PHI(l) 

PHI(2) 

g 

1 

90 

- 

O 

2 

0 

- 

0 

2 

90 

- 

2 

45 

- 

A 

3 

45 

-90 

3 

-45 

-90 

cQ> 

3 

45 

90 

<> 

A 

0 

90 

n 

4 

60 

-60 

4 

0 

60 

Figure  15  Representative  Inlet  Orientation 
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SIDE  ATTACHED  2D  INLET 


Figure  16  Inlet  Coordinate  System  Defined 
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Top  External  Chamfer 


Definitions: 


Location  Number 

Definition 

1 

Top  lip  leading  edge 

2 

Top  lip  chamfer  end 

3 

Midbody  end  and  boattail  start 

4 

Boattail  end  or  base 

5 

Bottom  lip  leading  edge 

6 

Forebody  maximum  height 

7 

Forebody  end  and  Midbody  start 

8 

Midbody  end  and  boattail  start 

9 

Internal  co^Rressitan  ramp  start 

10 

Internal  con— ii  tan  I'onp  end 

Figure  17  2-0  Inlet  Nomenclature  and  Assumptions 


NOMENCLATURE 


‘INLET 


lea6ing 

EDGE 

i- 


MAX'IMUM 

WIDTH 


.  vmimtrm 

I  I  ! 


BOATTAIL 

FAIRING 

START 


Vrrmntmttn 


*  ^ 

I - J-J 

-DIVERTER 


INPUTS 


ASSUMPTIONS 


Figure  19  Side  Attached  Inlet  Diverter  Nomenclature 
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FRONT  VIEW 


SIDE  VIEW 


— Z 


10 


10. 


Figure  20  Top  Attached  2-0  Inlet  Nomenclature 


'  >  BOATTAIL 

LEADING  MAXIMUM  FAIRING 

EDGE  WIDTH  START 

t  4-  + 


Figure  21  Top  Attached  Inlet  Diverter  Nomenclature 
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DEFINITIONS: 


LOCATION  NUMBER 

DEFINITION 

1 

CENTERBODY  NOSE  TIP 

2 

CENTERBODY  INTERMEDIATE  LOCATION 

3 

CENTERBODY  END 

4 

LIP  LEADING  EDGE 

5 

LIP  INTERMEDIATE  LOCATION 

6 

LIP  END  AND  MIDBODY  START 

7 

MIDBODY  END  AND  BOAHAIL  START 

8 

BOATTAIL  END  OR  BASE 

Figure  22  Axlsymmetric  Inlet  Nomenclature 
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ASSUMPTIONS 


Figure  23  Axisyinmetric  Inlet  Inputs 


Figure  24  Axisymmetric  Diverter  Nomenclature 
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NAMELIST  ARBBOD 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

XTIP 

- 

The  distance  along  the  x-axis  from 
the  origin  to  the  tip  of  the  forebody 

DTIP 

- 

The  width  of  the  forebody  at  the  tip 

ALENTH 

- 

The  length  of  the  forebody,  measured 
from  the  tip  to  the  base 

BLENTH 

- 

The  length  of  the  forebody,  measured 
from  the  base  of  the  nose  segment  to  the 
base  of  the  centerbody 

STACNT 

— 

The  number  of  stations  the  nose  segment 
will  be  divided  into  when  calculating 

Cm(x  up  to  a  maximum  of  nine 

SEGCNT 

“ 

The  number  of  linear  and  elliptic  arc 
segments  used  to  descibe  the  right 
half  of  the  forebody  base  cross  sectbn 

TTYPE 

9 

0  indicates  the  corresponding 
surface  segment  will  converge  to  a  point 
at  the  tip.  1  indicates  the  corresponding 
surfacebase  segment  will  converge  to  a 
line  segment  at  the  tip. 

BTYPE 

9 

1  -  point  (unused) 

2  ■  line  segment 

3  •  elliptic  arc  segment 

YSTART, 

ZSTART 

9 

The  coordinates  of  the  starting  point 
of  a  segment  in  the  base  cross  section 

YFINAL, 

ZFINAL 

9 

The  coordinates  of  the  end  point  of  a 
segment  in  the  base  cross  section 

YCENTR, 

ZCENTR 

9 

The  coordinates  of  the  center  of  an 
elliptic  arc  segment 

YLENTH, 

ZLENTH 

9 

The  semimajor  and  semiminor  axis  lengths. 

Note:  YLENGTH  does  not  have  to  be  the 
semimajor  axis.  For  a  circle,  only  YLENTH 
needs  to  be  specified. 

DSTART 

9 

Angle  of  the  starting  point  measured  from 
the  center  of  the  starting  ellipse, 
relative  to  the  positive  y-axis. 

□FINAL 

9 

Angle  of  the  end  point  measured  from  the 
center  of  the  ellipse  to  the  end  point  relative 
to  the  positive  y-axis. 

Figure  25  Arbitrary  Body  Inputs 
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CASE 10  ARBITRARY  BODY 

NAMELIST 

NO  1>T 

SFLTCON 


NALPHA-  2., 

ALPHA  -  0., 

1 . . 

NMACH  -  4. . 

MACH  -  2.50, 

3.00, 

REN  -  2.5E06. 

2.5E06, 

IREFQ  XCG-0. .  LREF 

Sarbboo 

-  3.333 

STACNT-  9.000, 

3. SC. 
2.5Ee6. 


4.ee, 

2.5Ee6. 


ALENTH-  1.111, 
BLENTH-  2.222. 
OTip  -  e.eee, 
BTYPE  -  3.. 

'TYPE  -  0.. 
osTART*-9e.ee0e, 
OFINAL-  00.0000. 
YCENTR-  0.0000, 
ZCENTR-  0.0000, 
YLENT1*-  0.1667, 
NEXT  CASE 


Figure  26  Cone/Cylinder  Input  Case 
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CASEIO  ARBITRARY  BODY 

NAMELIST 

NO  LAT 

IFLTCON 

NALPHA-  2.. 


ALPHA  -  e. . 

NMACH  -  4.. 

1.. 

MACH  -  2.5e, 

3.00. 

3.50, 

4.00. 

REN  -  2.5Ee6. 

2.5Ee6. 

2.5Ee6, 

2.5Ee6, 

SREFQ  XCG-e..  LREF  -  3.333. 
$ARB800 

STACNT-  9.eoe, 

SECCNT-  4.eee. 

ALENTH-  1.111, 

BLENTH-  2.222. 

OTiP  -  e.eee. 

$ 

BTYPE  -  2.. 

3.. 

2.. 

3. . 

TTYPE  -  0 . . 
YSTART-  0.0000. 
ZSTART— 0.1700, 

0.. 

0.. 

0., 

YFINAL-  0.1200, 

l.eE-30. 

0.0S78. 

1.0E-3e. 

ZF  INAL— 0.1700, 

1.0E-30. 

0.1367, 

1 .0E-30. 

YCENTR-  1.0E-30, 

0.1200, 

1 .0E-30, 

V  •  Ww  • 

ZCENTR-  1.0E-30, 

-0.1033. 

1 .0E-3e, 

0.1033, 

YLENT14-  1.0E-30. 

0.0667. 

1 .eE-30, 

0.0667, 

OF INAL-  l.eE-30. 
NEXT  CASE 

XA  AAAA 

1.0E-30. 

90.0000, 

Figure  27  Rounded  Triangular  Cone/Cylinder  Input  Case 


CASE ID  arbitrary  BODY 

NAMELIST 

NO  LAT 

lELTCON 


NALPHA-  2. . 

ALPHA  -  0. , 

NMACH  -  4. . 

1.. 

MACH  -  2.50, 

3.00. 

3.50, 

REN  -  2.5E06. 

2.5E06, 

2.5E06, 

SREFO  XCC-e. .  LREP 
lARBBOO 

STACNT-  9.000, 
SECCNT-  3.000, 
ALENTH-  .646. 

OTIP  m  0.1667. 

-  3.333, 

$ 

BTYPE  -  2., 

2.. 

2. , 

TTYPE  -  1 . . 
YSTART*  0.0000, 
ZSTART— 0.0000, 

0.. 

1., 

YPINAL-  0.2008, 

0.0833, 

0.0008, 

ZfINAL-  0.0000, 
NEXT  CASE 

0.0565, 

0.0565, 

4. Be, 

z.SEee,  $ 


Figure  28  Rhombus/Spatular  Nose  Input  Case 


CASEID  MISSILE  DATCOM  ARBITRARY  BODY  EXAMPLE 
DIM  FT 

SFLTCON  NMACH-2.e.  MACH-2.5.4.0.  REN-2.5E6.2.5E6. 

NALPHA-5.,  ALPHA-0.. 5.. 10.. 15.. 20.,$ 
IREFQ  XCG-1.667.$ 

$ARBBOO  XTIP-0.0,  OTIP-0.0. 

ALENTeVI .  1 1 1 ,  BLENT1^2.222. 

STACNT-9.0,  SEGCNT-2.0. 

TTYPE-0.0,  0.0, 

BTYPE-2.0.  3.0, 

YCENTR^1.E-30,  0.0, 

2CENTR-1 .E-30,  0.0, 

YLENTF^I .E-30,  .2886, 

2LENTF^1.E-30.  .1924, 

YSTART-0.0,  1.E-30. 

ZSTART-0.0.  1.E-30. 

YF INAL-. 2886.  1.E-30. 

ZFINAL-0.0.  1.E-30. 

OFINAL-90.0,$ 

SOSE 

NEXT  CASE 


Figure  29  Flat  Bottom  Ellipse/Centerbody  Input  Case 
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NAMELIST  EXPR 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

r  "  "  ■  " 

DEFINITION 

MACH 

1 

MACH  NUMBER 

NALPHA 

1 

NO.  ANGLES  OF  ATTACK  (2-20) 

ALPHA 

20 

ANGLES  OF  ATTACK  FOR  DATA 

SREF  > 

1 

REFERENCE  AREA  FOR  DATA 

LREF 

1 

LONGITUDINAL  REF.  LENGTH  FOR  DATA 

LATREF 

1 

LATERAL  REF.  LENGTH  FOR  DATA 

XCG 

1 

LONGITUDINAL  C.G.  FOR  DATA 

ZCG  ) 

1 

VERTICAL  C.G.  FOR  DATA 

CONF 

1 

CONFIGURATION  FOR  DATA 

SELECT  ONE  OF  THE  FOLLOWING 

BODY  -  BODY 

F1  -  WING 

F2  -  TAIL 

F3  -  3RD  FIN  SET 

F4-4THFIN  SET 

BF1  -  BODY-WING 

BF12- BODY-3  FIN  SETS 

BF1 234  -  BODY-4  FIN  SETS 

CN 

20 

C^  DATA  VS  ALPHA 

CM 

20 

C^,  DATA  VS  ALPHA 

CA 

20 

Ca  data  vs  ALPHA 

CY 

20 

Cy  DATA  VS  ALPHA 

OSN 

20 

Cn  DATA  VS  ALPHA 

CSL 

20 

1 

Ci  DATA  VS  ALPHA 

0  FOR  DEFAULT  VALUES,  SEE  NAMELIST  REFQ. 

Figure  30  Experimental  Data  Inputs 


Table  5  Airfoil  Designation  Using  the  NACA  Control  Card 

Ml  1-1 4099 


INPUT  NACA 
DESIGNATION 

NACA  SERIES 
AIRFOIL 

RESTRICTIONS 

0012.25 

4-Digit 

None.  Fractional  thickness  may  be  specified. 

23118.50 

5 -dig  it 

None.  Fractional  thickness  may  be  specified. 

2406-32 

4- Digit 
modified 

Sixth  digit  specifies  position  of  maximum 
thickness,  (%chord/l0),  and  must  be  a 

2,  3,  4,  5,  or  6. 

43006-65 

5-Digit 

modified 

Seventh  digit  specifies  postion  of  maximum, 

(%chord/10),  and  must  be  a  2, 3, 4, 5,  or  6. 

16-212.25 

1 -Series 

Second  digit  specifies  location  of  minimum  pressure, 
(%chord/10),  and  must  be  a  6, 8,  or  9.  Fractional 
thickness  may  be  specified. 

64-005 

64-205  A  =  0.6 
63A005 

652A215  A  =  0.8 
65,2A215  A  =  0.8 

6-Series 

Second  digit  specifies  location  of  minimum  pressure, 
(%chord/10),  and  must  be  a  3, 4, 5  or  6.  The  mean 
line  parameter  (A  =  xx)  must  be  a  decimal  between 

0.1  and  1 .0  (Defalut  is  1 .0). 

See  Note  1 . 

3-30.02.5-40.1 

A  B  C  D 

Supersonic 

See  Note  2. 

A  -  Section  type:1  =  Double  Wedge 

2  =  Circular  Arc 

3  =  Hexagonal 

8  -  Distance  from  leading  edge  to  position  of 
maxiOKim  thickness,  %  of  chord. 

C  -  Maximum  thickness,  %  of  chord. 

D  -  For  hexagonal  sections,  length  of  surface  of 
constant  thickness,  %  of  chord. 

Note  1  ■  The  program  does  not  distinguish  between  a  64,  2-220  and  a  64-220  specification. 
The  difference  in  coordinates  between  the  two  is  negligible. 


Note  2.  All  parameters  can  be  expressed  to  0.1%.  The  delimeter must  be  used. 
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Figure  31  Typical  "Stacked"  Case  Setup 
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Figure  32  "Configuration  Incrementing"  Case  Setup 
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4.0  OUTPUT  DESCRIPTION 


This  section  describes  the  types  of  output  available  from  the  code.  In 
most  cases  the  available  output  is  user  selectable,  that  is,  it  is  not  nor¬ 
mally  provided  and  must  be  specifically  requested  using  a  specialized  control 
card.  This  feature  permits  the  user  to  tailor  the  code  output  to  fit  his  par¬ 
ticular  application  without  extensive  reprogramming.  This  allows  him  to  find 
the  output  that  he  is  interested  in  without  having  to  wade  through  output  that 
does  not  interest  him. 

The  following  four  types  of  output  are  available  from  the  code: 

0  Nominal  output  -  This  output  is  always  provided  by  the  code  and  con¬ 
sists  of  output  from  the  input  error  checking  module  (CONERR),  a 
listing  of  the  inputs  for  each  case,  and  the  final  aerodynamic  re¬ 
sults  for  the  configuration. 

0  Partial  output  -  This  output  details  the  configuration  geometry  and 
the  intermediate  aerodynamic  calculations.  Special  control  cards  are 
available  so  that  the  user  can  select  the  quantity  and  types  of  out¬ 
put  desired. 

0  External  data  files  -  This  output  permits  the  user  to  create  external 
data  files  which  can  the  used  in  post-processing  programs,  such  as 
plotting  or  trajectory  programs.  Both  fixed  and  user  defined  format 
data  files  can  be  created  with  the  addition  of  simple  control  cards. 

0  Array  dumps  and  extrapolation  messages  -  This  output  permits  the  user 
to  print  internal  data  arrays  (DUMP)  or  to  investigate  design  chart 
extrapolations  during  the  course  of  the  case  execution  (PRINT  EXTRAP). 

A  special  feature  of  the  code  is  the  ability  to  utilize  it  as  a  subroutine 
package  to  a  larger  missile  design  code.  This  feature  is  described  in  Section 
4.5. 
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The  remainder  of  the  section  describes  each  of  these  output  data.  Exam¬ 
ples  of  each  output  page  are  also  Included  and  were  treated  from  the  example 
problems,  described  In  the  Appendix,  which  can  be  used  as  a  model  for  setting¬ 
up  another,  similar  configuration  or  be  used  as  a  means  to  check  the  proper 
operation  of  the  code. 

4.1  Nominal  Output 

Without  the  use  of  any  program  options  the  code  will  provide  three  types 
of  output.  First,  an  analysis  by  the  input  error  checking  routine  is  provid¬ 
ed.  It  lists  all  input  cards  provided  by  the  user  and  identifies  any  input 
errors  detected.  Second,  a  listing  of  all  input  cards,  grouped  by  case,  are 
provided;  included  in  this  output  is  an  error  analysis  from  the  major  input 
error  routine  MAJERR.  Finally,  the  total  configuration  aerodynamics  are  pro¬ 
vided  in  summary  form;  one  page  of  aerodynamic  output  is  supplied  for  each 
Mach  number  specified.  The  MAJERR  results  and  the  total  configuration  aero¬ 
dynamics  results  are  listed  in  succession  for  each  case, 

4.1.1  Input  Error  Checking 

The  purpose  of  the  input  error  checking  module  is  to  provide  single  pass 
error  checking  of  all  inputs.  If  an  error  is  detected,  it  is  identified  and 
an  appropriate  error  message  provided.  The  error  messages  are  designed  to  be 
self-explanatory.  In  some  cases,  errors  are  automatically  corrected  by  the 
routine,  although  the  routine  was  not  designed  to  be  a  comprehensive  error 
correction  utility. 

The  following  errors  are  automatically  corrected  by  the  code: 

0  No  terminating  comma  on  a  namelist  input  card 

0  No  terminating  or  "SEND"  on  a  namelist  input  ("8i"  on  IBM  systems) 

0  No  terminating  NEXT  CASE  for  the  case  inputs  for  single  case  or  last 
case  inputs. 
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Errors  detected  by  the  error  checking  routine  are  considered  either 
"FATAL"  or  "NON-FATAL",  a  "FATAL"  error  is  one  which  will  cause  the  code  to 
terminate  execution  abnormally;  examples  of  "FATAL"  errors  include  incorrect 
spelling  of  any  namelist  name,  incorrect  spelling  of  any  variable  name,  and 
any  drastic  input  error  in  a  namelist  input,  such  as  leaving  out  an  equals 
sign  in  a  constant  definition.  All  "FATAL"  errors  are  clearly  identified  on 
the  output.  A  "NON-FATAL"  error  is  one  which  will  not  cause  the  program  to 
terminate  execution;  an  example  of  a  "NON-FATAL"  error  is  leaving  off  the  de¬ 
cimal  point  on  numeric  constants  all  Missile  Datcom  inputs  are  either  REAL  or 
LOGICAL  regardless  of  the  variable  name  assigned.  "NON-FATAL"  errors  will  not 
cause  the  code  to  stop  execution,  whereas,  "FATAL"  errors  will  cause  the  code 
to  stop  execution  after  input  error  checking  has  been  completed. 

An  example  output  from  CONERR  is  shown  in  Figure  33.  This  figure  illus¬ 
trates  the  array  of  input  errors  checked  by  CONERR.  Several  additional  fea¬ 
tures  of  the  output  are  as  follows; 

0  All  user  defined  input  cards  are  assigned  a  sequential  "line  num¬ 
ber".  This  serves  as  means  of  inputs  from  the  code  generated  ones 

(all  code-created  input  cards  are  not  identified  with  a  "line  num¬ 
ber").  This  scheme  also  permits  the  user  to  quickly  identify  those 
input  cards  in  error  so  that  efficient  correction  of  input  errors  can 
be  performed. 

0  All  input  cards  are  listed  as  input  by  the  user.  To  the  right  of 
each  input  card  is  a  listing  of  any  errors  encountered  in  processing 
that  card.  If  no  such  error  message  appears  then  the  input  was  in¬ 
terpreted  as  being  correct. 

0  In  many  cases  alphanumeric  constants  are  available  (see  Table  3). 
Hence  the  user  does  not  need  to  memorize  a  numeric  scheme  of 
"flags".  Since  some  computers  do  not  recognize  alphanumeric  con¬ 
stants  as  namelist  constants,  they  are  automatically  converted  by  the 
code  to  their  numeric  equivalent.  A  message  is  printed  to  identify 
the  substitutions  performed.  The  example  input  in  Figure  34  shows 
replacements  for  CONE  and  OGIVE. 
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In  order  to  permit  column  Independent  Inputs  the  code  will  automatically 
adjust  some  of  the  input  cards  to  begin  in  columns  1  or  2.  All  control  cards 
will  be  automatically  shifted  to  start  in  column  1;  all  namelists  which  begin 
in  column  1  will  be  shifted  to  column  2.  If  any  input  card  cannot  be  shifted 
to  conform  to  this  scheme,  an  error  message  will  be  produced.  As  a  general 
rule,  column  80  of  namelist  inputs  should  be  left  blank  so  that  the  code  can 
shift  the  card  image,  if  necessary, 

4.1.2  Listing  of  Case  Input  Data 

Figure  34  shows  the  first  page  of  outputs  for  a  case  without  CONERR  de¬ 
tected  errors.  Then  Figure  35  shows  the  next  page  of  output  which  lists  all 
input  cards  for  the  case  (down  to  the  NEXT  CASE  control  card).  If  the  input 
for  a  case  is  from  a  previous  case  (through  use  of  the  SAVE  control  card)  only 
the  new  case  inputs  are  listed.  All  saved  inputs  are  not  repeated  in  subse¬ 
quent  case  input  summaries. 

After  the  case  data  have  been  read,  the  data  setup  for  the  case  is  ana¬ 
lyzed  by  the  case  major  error  checking  module  (MAJERR).  The  purpose  of  this 
second  error  checking  is  to  insure  that  the  data  input,  although  syntax  error 
free,  properly  defines  a  case  to  be  run.  Examples  of  errors  detected  in 
MAJERR  include  valid  flight  condition  inputs,  valid  reference  condition  in¬ 
puts,  and  that  geometry  has  been  defined.  In  most  cases  errors  detected  by 
MAJERR  are  corrected  with  assumed  defaults.  If  any  MAJERR  error  message  is 
produced,  the  user  should  verify  the  "fix-up"  taken  by  the  code.  In  some 
cases  a  "fix-up"  is  not  possible;  an  appropriate  error  message  and  a  sugges¬ 
tion  for  correcting  the  error  is  provided.  If  a  "fix-up"  is  not  possible  the 
case  will  not  run. 

4.1.3  Case  Total  Configuration  Aerodynamic  Output  Summary 

As  shown  in  Figure  36,  the  total  configuration  aerodynamics  are  provided 
in  compact  form  for  easy  review.  The  aerodynamics  are  summarized  as  a  func¬ 
tion  of  angle  of  attack  (ALPHA)  in  the  user  specified  system  of  units,  the 
nomenclature  is  as  follows: 
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CN  -  Normal  force  coefficient 

CM  -  Pitching  moment  coefficient 

CA  -  Axial  force  coefficient 

CY  -  Side  force  coefficient 

CLN  -  Yawing  moment  coefficient 

CLL  -  Rolling  moment  coefficient 

CNA  -  Normal  force  coefficient  derivative  with  ALPHA 

CMA  -  Pitching  moment  coefficient  derivative  with  ALPHA 

CYB  -  Side  force  coefficient  derivative  with  BETA 

CLNB  -  Yawing  moment  coefficient  derivative  with  BETA 

CLLB  -  Rolling  moment  coefficient  derivative  with  BETA 

CL  -  Lift  coefficient 

CD  -  Drag  coefficient 

CL/CD  -  Lift  to  drag  ratio 

XCP  -  Center  of  pressure  from  the  moment  reference  center  divided  by 

reference  length 
DX/L  -  (XCP  -  XCG)/LREF 

All  coefficients  are  based  upon  the  reference  areas  and  lengths  specified 
at  the  top  of  the  output  page.  The  derivatives  CNA  and  CMA  are  computed  by 
numeric  differentiation  of  the  CN  and  CM  curves,  respectively;  precise  deriva¬ 
tives  are  only  obtained  when  the  angle  of  attack  range  specified  is  narrow. 
The  derivatives  CYB,  CLNB  and  CLLB  are  determined  by  perturbing  the  sideslip 
angle  by  one  degree,  recalculating  the  configuration  forces  and  moments,  and 
then  differencing  with  the  user  specified  orientation.  Hence,  the  longitud¬ 
inal  and  lateral  derivatives  will  probably  not  be  numerically  identical  for 
those  conditions  which  should  produce  identical  results  if  they  were  both  cal¬ 
culated  by  the  same  method.  ^ 

A  significant  decrease  in  computational  time  is  realized  when  the  calcula¬ 
tion  of  lateral-directional  derivatives  are  suppressed  using  the  control  card 
NO  LAT.  For  these  cases,  the  CYB,  CLNB,  and  CLLB  data  fields  are  filled  with 
blanks. 
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When  selecting  TRIM,  the  output  is  provided  in  a  form  similar  to  Figure 
37.  When  running  a  trim  case  the  derivatives  due  to  ALPHA  and  BETA  are  not 
available.  The  panels  which  were  deflected  to  trim  the  configuration  are  in¬ 
dicated  by  the  "VARIED"  citation  next  to  them. 

The  format  for  the  values  of  the  numbers  in  the  printed  output  has  been 
assumed  based  on  typical  magnitudes  for  missile  aerodynamic  coefficients.  In 
some  cases,  a  user  specified  reference  area  and/or  length  will  cause  the  re¬ 
sults  to  underflow  or  overflow  the  formal  selected.  For  these  cases  the  user 
should  adjust  his  reference  quantities  by  powers  of  ten  to  get  the  data  to  fit 
the  format  specified. 

4.2  Partial  Output 

Partial  output  consists  of  geometry  calculation  details,  intermediate 
aerodynamic  results,  or  auxiliary  data,  such  as  pressure  distributions.  Each 
of  these  output  types  are  printed  through  the  addition  of  control  cards  input 
for  each  case.  In  all  cases,  partial  output  requested  for  one  case  is  not 
automatically  selected  for  subsequent  cases,  and  the  control  cards  must  be 
re-input.  This  permits  the  user  to  be  selective  on  the  amount  and  types  of 
output  desired. 

A  special  control  card  PART  permits  the  user  to  request  all  geometric  and 
aerodynamic  partial  output.  Due  to  the  amount  of  output  produced,  this  option 
should  be  used  sparingly  or  when  details  of  the  calculations  are  desired. 

There  is  one  geometry  for  which  no  partial  output  is  obtained.  This  is 
when  bodies  with  arbitrary  cross  sections  are  input. 

The  following  paragraphs  describe  the  output  received  when  partial  output 
is  requested. 

4.2.1  Geometric  Partial  Output 

Details  of  the  geometry  are  provided  when  the  PART  or  PRINT  GEOM  control 
cards  are  included  in  the  case  inputs.  Figure  38  shows  the  output  created 
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when  the  PRINT  6E0M  BODY  control  card  is  used.  Detailed  are  the  results  of 
the  geometric  calculations  for  the  body.  Included  are  such  items  as  planform 
area,  surface  (wetted)  area,  and  the  mold  line  contour. 

If  fins  are  present  on  the  configuration,  two  types  of  fin  geometry  data 
are  produced  when  PRINT  GEOM  FINl  or  PART  is  requested.  As  shown  in  Figure 
39,  the  description  of  the  panel  airfoil  section  is  provided.  Following  that, 
shown  in  Figure  40,  is  a  summary  of  the  major  geometric  characteristics  of 
such  planform;  note  that  fin  planform  geometry  data  is  given  for  one  panel  of 
each  fin  set,  since  it  is  assumed  that  each  fin  of  a  fin  set  is  identical.  If 
a  panel  is  made  up  of  multiple  segments,  the  geometric  data  is  provided  by 
panel  segment  (each  segment  is  assigned  a  number  starting  at  the  root).  Total 
panel  set  of  characteristics  is  also  provided.  This  total  panel  data  repre¬ 
sents  an  equivalent  straight-tapered  panel,  which  is  used  for  most  of  the 
aerodynamic  calculations.  The  thickness-to-chord  ratio  shown  for  each  segment 
is  that  value  at  the  segment  root;  for  the  total  panel,  it  is  an  "effective" 
value. 

If  an  airbreathing  inlet  is  specified  (presently  available  only  at  super¬ 
sonic  speeds)  the  output  is  similar  to  that  in  Figure  41.  This  output  re¬ 
flects  the  user  input  definition  for  the  inlet  design  specified.  It  is  pro¬ 
vided  if  the  PRINT  GEOM  INLET  or  PART  control  cards  are  included  in  the  input 
case. 

4.2.2  Aerodynamic  Partial  Output 

The  output  on  the  configuration  aerodynamics  is  most  extensive  when  PRINT 
AERO  or  PART  is  specified.  As  shown  in  Figures  42,  43,  44  and  45,  output  is 
created  for  the  body  and  each  fin  set  on  the  configuration.  In  addition,  for 
any  subsonic/transonic  Mach  number  (less  than  1.4)  an  analysis  by  the  Airfoil 
Section  Module  is  made,  which  involves  a  potential  low  analysis  of  the  airfoil 
section  using  conformal  mapping  is  made.  Figure  46  shows  the  output  produced 
from  this  analysis.  If  an  inletted  configuration  additional  partial  output  is 
included  to  summarize  the  inlet  external  aerodynamics.  Figure  48  presents  an 
example  of  this  output. 
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As  shown  1n  Figure  42,  the  body  alone  partial  aerodynamic  output  for  nor¬ 
mal  force  lists  the  axial  force  contributors,  potential  normal  force  (CN- 
POTENTIAL),  viscous  normal  forces  (CN-VISCOUS) ,  potential  pitching  moment  (CM- 
POTENTIAL),  viscous  pitching  moment  (CM-VISCOUS) ,  and  the  crossflow  drag  co¬ 
efficient  (CDC).  The  cross-flow  drag  proportionality  factor  at  subsonic  and 
transonic  speeds  is  also  given  for  reference.  These  data  are  similar  to  that 
obtained  for  elliptical  bodies. 

Figure  43  details  the  fin  normal  force  calculations  by  fin  set.  Each 
panel's  contribution  to  the  configuration  normal  force  is  described.  The 
column  titled  CN-POTENTIAL  is  the  potential  contribution  and  the  column  titled 
CN-VISCOUS  is  the  viscous  contribution.  Their  sum  is  given  in  the  column 
titled  CN-TOTAL.  CNAA  is  the  nonlinear  variation  of  normal  force  due  to  angle 
of  attack  and  ALPHA  EQUIV  is  the  panel  angle  of  attack  due  to  its  roll  posi¬ 
tion  on  the  body.  Figure  45  illustrates  the  fin  axial  force  contributors  and 
Figure  46  presents  an  example  of  the  fin  pitching  moment  contributors. 

The  analysis  by  the  Airfoil  Section  Module  is  provided  in  a  format  similar 
to  Figure  46.  If  any  mach  number  specified  produces  supersonic  flow  on  the 
airfoil  surface,  the  message  "CREST  CRITICAL  MACH  NUMBER  EXCEEDED"  will  be 
printed;  approximation  of  the  airfoil  section  data  is  then  assumed.  These  fin 
aerodynamic  increments  are  repeated  for  each  fin  set  on  the  configuration. 
Note  that  the  Airfoil  Section  Module  assumes  that  the  panels  have  sharp  trail¬ 
ing  edges.  Any  panel  input  with  a  non-sharp  trailing  edge  will  have  its  aero¬ 
dynamic  characteristics  set  as  though  the  airfoil  was  "ideal".  This  assump¬ 
tion  is  approximate  for  preliminary  design. 

Figure  47  shows  the  aerodynamic  output  available  when  inlets  are  specified 
on  the  configuration.  It  is  provided  when  PRINT  AERO  INLET  or  PART  is  speci¬ 
fied  in  the  case  Inputs.  The  aerodynamics  summarized  for  inlets  assume  that 
the  inlets  are  operating  at  their  critical  Mach  number.  No  internal  aerody¬ 
namics  are  calculated  or  included  in  the  output. 

After  the  aerodynamic  details  for  each  component  of  the  configuration  are 
output,  the  aerodynamic  calculations  for  the  synthesis  of  the  complete  config¬ 
uration  follows.  As  shown  in  Figure  48,  the  aerodynamics  for  the  body  are 
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first  presented.  This  is  followed  by  the  following  output  for  each  fin  set. 
In  the  case  of  the  example  case  ^1n  set  1  results  would  be  followed  by  fin  set 
2  results  for  each  of  the  following  outputs: 


0  "FIN  SET  PRESENCE  OF  THE  BODY"  -  This  summarizes  the  aerodynamic  in- 
crementals  of  the  most  forward  set  of  fins  with  the  influence  of  the 
body.  Figure  49  presents  the  example  of  this  output.  The  left-most 
six  columns  include  the  effect  of  body-on-fIn  component  interfer¬ 
ence.  The  right-most  columns  represent  the  contribution  to  each 
panel  to  configuration  aerodynamics,  and  include  the  effect  of  body- 
on-fin  interference,  these  values  are,  in  effect,  individual  panel 
louis.  The  panel  characteristic  values  included  are  "AEQn"  (the 

panel  equivalent  (local)  angle  of  attack)  and  "CNn"  (the  panel  normal 

force  coefficient).  The  sign  convention  is  as  follows:  a  positive 

panel  normal  force,  hence,  equivalent  angle  of  attack,  produces  a 
negative  roll  moment.  Hence,  panels  on  the  right  side  of  the  con¬ 
figuration  will  produce  loads  and  angles  of  attack  opposite  in  sign 

to  those  on  the  left  side  of  the  configuration  even  though  they  pro¬ 
duce  the  same  physical  force  loading. 

0  "BODY-FIN  SET"  -  Aerodynamics  for  the  body  plus  most  forward  set  of 
fins  configuration.  It  is  produced  through  addition  of  the  body 
alone  and  wing  in  presence  of  the  body  incrementals,  described 
above.  The  results  shown  in  Figure  50,  Include  the  component  carry¬ 
over  factors  K-W(B)  (wing  in  presence  of  the  body  carryover  due  to 
angle  of  attack),  K-B(W)  (body  in  presence  of  the  wing  carryover  due 
to  angle  of  attack),  KK-B(W)  (body  in  presence  of  the  wing  carryover 
due  to  panel  deflection),  JCP-W(B)  (wing  in  presence  of  the  body 
carryover  center  of  pressure),  and  XCP-B(W)  (body  in  presence  of  the 
wing  carryover  center  of  pressure).  This  output  is  repeated  for  the 
body  plus  the  first  and  second  most  forward  lifting  surfaces,  if  one 
exists  on  the  configuration.  This  example  includes  two  fin  sets  so 
the  next  page  of  partial  output  would  look  like  Figure  51.  If  addi¬ 
tional  fin  sets  are  present  on  the  configuration  additional  pages  are 
output  with  each  one  successively  included. 
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"CARRYOVER  INTERFERENCE  FACTORS"  -  This  page  of  partial  output  sum¬ 
marizes  the  carryover  factors  listed  in  the  paragraph  above.  These 
were  included  in  the  body  plus  fin  set  calculations.  An  example  of 
this  output  is  presented  in  Figure  52. 

0  "COMPLETE  CONFIGURATION"  -  Complete  configuration  aerodynamics.  This 
output  is  illustrated  in  Figure  36.  The  values  are  obtained  by  sum¬ 
ming  the  body-wing  and  tail  in  the  presence  of  the  wing  flow  field 
data. 

If  the  PRINT  AERO  BEND  or  PART  control  card  is  used,  the  code  will  compute 
and  print  panel  bending  moment  coefficients  for  each  fin  set  on  a  separate 
page.  One  page  is  shown  in  Figure  53.  The  sign  convention  is  that  assumed 
for  the  individual  panel  loads  and  equivalent  angles  of  attack,  noted  above. 
The  bending  moment  coefficients  are  based  upon  the  reference  area  and  longi¬ 
tudinal  length  given  at  the  top  of  the  page.  The  moments  are  referenced  about 

the  fin-body  structure  specified  by  the  root  chord  span  station. 

Figure  54  illustrates  the  panel  hinge  moments  coefficients  computed  when 
the  control  cards  PRINT  AERO  HINGE  or  PART  are  used.  The  reference  area  and 

longitudinal  reference  length  given  at  the  top  of  the  page  are  used.  All  mo¬ 

ments  are  computed  about  the  hinge  line,  which  is  defined  using  namelist 
DEFECT. 

If  TRIM  is  specified,  the  user  can  selectively  print  the  six  untrimmed 

static  aerodynamic  tables  used  in  the  trim  process.  An  example  is  shown  in 

Figure  55.  The  code  computes  the  six-component  aerodynamics  at  ten  deflection 

angles  for  each  specified  angle  of  attack,  then  interpolates  for  C  =0.  Note 

m 

that  this  trim  process  can  be  used  to  create  control  authority  data,  effec¬ 
tively  giving  the  user  10  deflection  angles,  20  angles  of  attack,  and  20  Mach 
numbers  per  input  case. 

4.2.3  Pressure  Distribution  Data 

If  the  Mach  number  is  supersonic  (M>1.2),  the  user  has  the  option  to  print 
the  surface  pressure  distributions  over  the  body  and  fins.  This  option  is 
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selected  only  through  the  addition  of  the  control  card  PReSSURES.  Since  three 
body  alone  supersonic  methods  are  available  (Van  Dyke  Hybrid,  Second-Order 
Shock  Expansion  (SOSE),  and  Newtonian  flow)  the  capability  exists  to  output 
the  pressure  distribution  data  from  any  one  of  these  methods.  The  method  to 
be  used  in  the  calculation  of  the  pressure  data  is  controlled  with  the  control 
cards  SOSE  and  HYPER;  if  neither  control  card  is  input,  the  Van  Dyke  Hybrid 
method  is  selected.  Because  of  the  nature  of  the  calculations,  body  alone 
pressure  are  printed  for  angles  of  attack  less  than  or  equal  to  15  degrees 
when  using  the  Hybrid  or  SOSE  techniques. 

The  capability  also  exists  for  the  user  to  output  the  pressure  distribu¬ 
tion  data  over  fins  at  any  Mach  number  greater  than  1.05.  This  option  is  also 
controlled  by  the  PRESSURES  control  card.  Due  to  the  nature  of  the  method, 
only  pressure  distribution  data  at  zero  angle  of  attack  is  presently  output. 

Figures  56,  51  and  58  illustrate  typical  output  produced  when  PRESSURES  is 
specified.  The  format  of  Figure  56  is  only  available  when  SOSE  is  specified; 
all  other  body  alone  pressure  methods  produce  output  similar  to  Figure  57  for 
bodies.  Figure  58  is  representative  of  fin  pressure  distribution  output. 
Note  that  calculation  of  pressures  is  a  time-consuming  process;  much  higher 
computational  times  will  be  required. 

All  body  pressure  distribution  data  is  based  on  a  configuration  that  has 
body  diameter  of  unity;  that  is,  the  configuration  is  expressed  in  calibers 
(or  body  diameters).  The  longitudinal  stations  at  which  pressure  coefficient 
data  is  desired  cannot  be  user  specified;  however,  sufficient  data  is  provided 
to  permit  accurate  interpolation  for  most  applications. 

4.3  Dynamic  Derivatives 

As  shown  in  Figure  59,  the  total  configuration  dynamic  derivatives  are 
provided  in  compact  form  for  easy  interpretation.  The  dynamic  derivatives  are 
summarized  as  a  function  of  angle  of  attack  in  the  user  specified  units.  The 
coefficients  provided  are  as  follows: 
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CNQ  -  Normal  force  coefficient  due  to  pitch  rate 

CNAD  -  Normal  force  coefficient  due  to  acceleration  in  angle  of 

attack  (a) 

CMO  -  Pitching  moment  coefficient  due  to  pitch  rate 

CHAO  -  Pitching  moment  coefficient  due  to  acceleration  in  angle  of 

attack  (a) 

Note:  For  body  alone  and  body  +  fin  set  data  CMQ  and  CHAD  are  presented  as 
the  sum  CMO+CMAD. 

The  dynamic  derivatives  are  printed  after  all  static  coefficients  and  par¬ 
tial  static  aerodynamics  are  printed.  If  a  build  or  part  card  is  input,  addi¬ 
tional  dynamic  derivatives  for  partial  configurations  and/or  configuration 
components  are  printed, 

4.4  External  Data  Files 


The  code  has  the  capability  to  be  used  in  conjunction  with  other  missile 
design  tools,  such  as  post-processing  plotting  programs  or  trajectory  pro¬ 
grams,  Fixed  format  aerodynamic  data  is  output  as  an  external  data  file  with 
the  addition  of  the  PLOT  control  card.  Included  in  this  data  file  are  the  six 
component  forces  and  moments  based  upon  the  user  specified  reference  quanti¬ 
ties.  In  order  to  print  component  buildup  data  to  the  plot  file  the  BUILD  and 
PLOT  control  cards  must  be  present  in  the  case.  Refer  to  Volume  2  for  a  com¬ 
plete  description  of  the  plot  file  output  format. 

An  option  to  create  a  user  specified  format  data  file  is  also  available. 
The  control  cards  WRITE  and  FORMAT  have  been  designed  for  easy  access  to  this 
capability. 
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4.5  Extrapolation  Messages  and  Array  Dumps 


As  shown  in  Figure  60,  the  extrapolation  messages  are  summarized  for  all 
design  charts  which  have  been  extrapolated  during  the  execution  of  the  case. 
Since  many  of  the  aerodynamic  methods  do  not  include  design  charts,  but  are 
either  closed-form  equations  or  complete  theoretical  methods,  this  option  is 
most  useful  in  the  subsonic  and  transonic  Mach  regimes.  Extrapolation  mes¬ 
sages  are  only  provided  if  the  control  card  PRINTE  EXTRAP  appears  in  the  case 
inputs.  The  data  titled  "ROUTINE  TRACE-BACK"  lists  the  dects  called  when  the 
look-up  was  performed;  "X"  represents  the  independent  variable  and  "Y"  repre¬ 
sents  the  dependent  variable  in  the  extrapolation. 

When  it  is  necessary  to  examine  the  values  stored  in  internal  data  arrays 
the  DUMP  control  card  can  be  used.  This  control  card  causes  the  contents  of 
the  named  data  arrays  to  be  printed  in  a  form  similar  to  Figure  61.  Array 
dumps  are  provided  for  each  Mach  number  of  the  input  case,  and  represent  the 
data  block  contents  at  aerodynamic  calculation  completion.  Volume  2  of  this 
User's  Guide  details  the  arrays  and  their  contents. 

Note  that  all  data  arrays  are  Initialized  to  a  constant  named  "UNUSED", 

-30 

which  is  preset  to  a  value  of  1  x  10  .  Hence,  any  array  element  which 
contains  this  constant  was  not  changed  during  execution  of  the  case  (since  it 
is  highly  unlikely  that  this  constant  will  result  from  any  calculation).  This 
scheme  permits  rapid  "tracking"  of  program  calculation  sequences  while  in  "de¬ 
bug"  mode. 

4.6  Special  Program  Capabilities 

The  code  is  currently  being  used  both  as  a  stand-alone  aerodynamic  predic¬ 
tion  code  and  as  a  subroutine  package  to  a  larger  missile  design  code.  Since 
all  of  the  printed  output  is  controlled  by  ^int  flags,  some  of  which  are  not 
directly  user  selectable,  the  code  can  be  used  as  part  of  a  larger  system  rel¬ 
atively  easily.  Users  desiring  to  utilize  this  capability  should  contact  the 
procuring  agency  for  instructions. 
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On  computer  systems  which  do  not  utilize  the  namelist  Input  scheme,  a 
fixed-format  Input  mode  can  be  easily  treated.  The  procuring  agency  can  pro¬ 
vide  advice  on  this  option  when  necessary. 
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Figure  34  Case  Input  Listing 
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Figure  35  Example  of  Default  Substitutions  for  Incomplete  Case  Inputs 
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Figure  36  Total  Configuration  Aerodynamic  Output  Summary 
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Figurt  37  Trlnmed  Output  Summary 


101 


Figure  38  Body  Geometry  Output 
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Figure  39  Airfoil  Geometry  Output 
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Figure  40  Fin  Geometry  Output 
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Figure  41  Inlet  Geometry  Output 
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Figure  42  Body  Alone  Aerodynamic  Partial  Output 
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Figure  43  Fin  Nornwl  Force  Partial  Output 
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Figurt  44  F1n  Axitl  Force  Partial  Output 
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Figure  45  Fin  Pitching  Moment  Partial  Output 
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Figure  46  Airfoil  Section  Aerodynamics  Partial  Output 
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Figure  47  Inlet  Aerodynamic  Partial  Output 
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Figure  48  Body  Alone  Synthesis  Partial  Output 
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Figure  49  Fin  Set  In  Presence  of  the  Body  Partial  Output 


Figure  50  Body  Plus  Fin  Set  Partial  Output 
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Figure  52  Carryover  Interference  Factors  Partial  Output 
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Figure  53  Panel  Bending  Moment  Partial  Output 
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Figure  58  Fin  Pressure  Distribution  Output 


Figure  59  Dynamic  Derivative  Output 
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Figure  60  Extrapolation  Massage  Output 
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Figure  61  Internal  Array  Dump  Output 
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A.  EXAMPLE  PROBLEMS 


This  appendix  presents  two  example  problems  for  use  in  verifying  the  pro¬ 
per  operation  of  the  computer  code,  or  for  use  as  a  model  for  the  setup  of 
inputs  of  similar  configurations.  The  first  example  is  a  simple  tangent  ogive 
nose-cylinder  circular  body.  It  has  a  planar  wing  (two  panels)  and  a  cruci¬ 
form  set  of  tails  orientated  in  the  "plus"  configuration.  There  are  two  input 
cases  for  this  problem.  The  first  demonstrates  the  output  generated  when  the 
PART  control  card  is  included.  The  second  demonstrates  the  output  created 
when  trim  has  been  requested;  note  that  trim  partial  output  has  been  requested 
with  the  PRINT  AERO  TRIM  control  card. 

The  second  problem  is  a  body-tail -inlet  configuration.  Three  Mach  numbers 
are  requested,  one  subsonic,  one  transonic,  and  one  supersonic.  Although  all 
three  Mach  numbers  could  be  run  in  a  single  case,  they  have  been  divided  into 
three  separate  cases  to  illustrate  the  “SAVE"  feature  as  well  as  selecting 
particular  output  for  illustration.  The  inlet  computations  are  done  only  for 
the  supersonic  Mach  number  case,  since  the  methods  are  only  valid  super¬ 
sonically. 
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A.l  Example  Problem  1 


The  first  example  problem  is  shown  in  Figure  A-1 .  It  is  comprised  of  a 
3-caliber  tangent  ogive  nose  attached  to  a  cylindrical  body;  a  triangular 
monoplane  set  of  wings;  and  a  cruciform  set  of  tails  orientated  in  the  "plus" 
position.  The  first  case  is  a  simple  angle  of  attack  sweep;  component  build¬ 
up  data  and  partial  output  are  requested.  The  second  case  is  a  trim  of  the 
configuration  using  the  two  horizontal  tail  surfaces. 
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A. 2  Example  Problem  2 


The  configuration  for  this  example  is  sketched  in  Figure  A-2.  The  figure 
is  a  modified  copy  of  the  wind  tunnel  model  drawings  from  NASA  Technical  Memo¬ 
randum  84557.  The  model  definition  in  Figure  A-2  is  representative  of  the 
detail  normally  found  on  design  drawings. 

This  example  has  been  divided  into  subsonic,  transonic,  and  supersonic 
cases.  Each  case  is  run  for  one  Mach  number.  Although  all  three  Mach  numbers 
could  have  been  run  in  one  case  they  were  run  separately  to  demonstrate  the 
SAVE  capability. 

This  example  provides  a  check  case  for  the  inlet  option.  It  can  be  used 
by  the  user  to  make  sure  that  he  understands  the  inputs. 
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Inlet-inlet  fairing  details 
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Figure  A-2  (cont.) 
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